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SECTION  I 
INTRODUCTION 


Solution  of  the  supersonic  boundary  layer  equations  requires  accurate  modeLs 
to  describe  the  turbulent  transport  of  energy  and  momentum.  The  transport 
coefficients  cannot  be  derived  from  first  principles,  but  must  be  deduced  from 
well-designed  experimental  studies.  Ideally,  it  is  preferred  to  measure  the 
turbulence  terms  directly,  e.g.,  using  the  hot  wire  anemometer,  although 
documentation  of  the  mean  flow  field  is  still  necessary  to  interpret  the  hot 
wire  data.  Even  in  the  absence  of  direct  information,  the  so-called  "inverse" 
or  "indirect"  method  can  be  applied  to  extract  the  transport  properties  from 
detailed  mean  flow  measurements  by  means  of  the  time  averaged  conservation 
equations.  Although  considerable  progress  in  this  direction  has  been  made 
for  the  zero  pressure  gradient,  flat  plate  boundary  layer,  there  is  still  a 
dearth  of  specific  information  concerning  flows  with  a finite  continuous 
variation  in  static  pressure.  This  problem  arises  from  the  fact  that  while 
numerous  experimental  studies  have  been  made  of  boundary  layers  with  pressure 
gradients  (see,  e.g.,  the  exhaustive  catalog  of  experimental  results  compiled 
by  Fernholz  and  Finley^)  few  have  been  sufficiently  complete  and  reliable  to 
successfully  apply  the  "indirect"  method.  In  fact,  for  continuous  pressure 
gradient  flows,  no  direct  hot  wire  measurements  of  transport  properties  are 
available  and  only  Sturek“  has  reported  on  turbulent  shear  stresses  calculated 
from  mean  flow  data. 

The  present  study  was  undertaken,  therefore,  to  systematically  study  the  in- 
fluence of  a continuous  adverse  pressure  gradient  on  a two-dimensional,  super- 
sonic, adiabatic  wall  boundary  layer.  The  adverse  pressure  gradient  was 
generated  using  a curved  ramp,  located  on  the  test  section  floor  of  the  FACC 
Mach  3 wind  tunnel,  and  designed  to  produce  an  isentropic  compression  with 
constant  dp/dx.  Two  ramps  were  designed  for  this  purpose,  one  corresponding 
to  a weak  and  the  other  to  a moderate  pressure  gradient.  UTiile  it  may  have 
been  more  desirable  to  design  the  ramps  for  constant  pressure  gradient 
parameter  5^  = (dp/dx),  and  thus  maintain  a constant  influence  of  the 


pressure  gradient  on  the  boundary  Laver,  die  local  values  ot  6,  and  t could 
not  be  predicted  prior  to  the  experiment.  On  the  other  hand,  it  is  shown  later 
that  the  boundary  layer  is  in  1 state  ol  local  equilibrium  and,  therefore,  is 
controlled  only  by  local  conditions,  i.e.,  c. , Me , a|<,  etc.  Thus,  maintaining 
the  pressure  gradient  constant  allows  the  overall  flow  field  to  lie  characterized 
by  i single  parameter  dp/dx  (or,  more  exactly,  by  where  '"/Tw  is  evaluated 

immediately  upstream  of  the  ramp)  while  providing  the  opportunity  to  study  the 
effect  ot  the  local  and  in  this  sense  the  experiment  yields  greater 

inf ormat ion . 

The  ultimate  goal  ot  this  program  is  the  direct  measurement  of  the  turbulent 
shear  stresses  using  the  x- array  hot  wire  anemometer  In  preparation  for  this 
Cask,  detailed  measurements  ot  the  mean  tlow  tield  over  the  curved  ramps  have 
aeen  carried  out.  This  report  describes  these  measurements  and  the  results  of 
the  data  analysis.  Details  ot  the  experiment  are  presented  in  Section  2,  while 
the  method  ot  malvsis  is  summarized  in  Section  3 and  the  results  are  discussed 
in  Section  *.  In  particular,  it  is  shown  that  the  mean  flow  profiles,  when 
suitably  transformed,  agree  with  the  universal  "wall-wake"  velocity  profile 
ind  that  the  pressure  gradient  parameter  ‘ ^ correlates  the  data  with  low  speed 
results.  In  addition,  the  distribution  across  the  boundary  Layer  of  turbulent 
shear  stress,  mixing  length,  and  eddv  viscosity  were  determined  using  the 
"indirect " method  modified  by  Sturek-  to  account  for  the  effects  of  longitudinal 
curvature  . 


SECTION  II 
THE  EXPERIMENT 


1*  WIND  TUNNEL 

The  experiment  was  carried  out  in  the  FACC  Mach  3 supersonic  wind  tunnel  (SWT). 

This  is  a continuous  flow  facility  with  a 7.87  cm  by  8.b4  cm  test  section 

located  40 . b cm  downstream  of  the  throat  section.  The  turbulent  boundary  layer 

at  Che  entrance  to  the  cest  section  is  fully  developed  and  is  approximately 

0.7  cm  thick.  All  tests  were  conducted  for  stagnation  conditions  ol  .973  x 
S 2 

10  N/M  and  317° K corresponding  to  a nominal  Moo  = 3 and  unit  Reynolds  number 
= 6.57  x lO^/meter. 

2.  RAMP  DESIGN 

Using  Method  of  Characteristics  theory,  and  ignoring  the  effect  of  the  boundary 
layer,  two  continuously  curved  ramps,  designated  Ramp  1 and  Ramp  3,  were  de- 
signed to  produce  two-dimensional,  constant  pressure  gradient  (dp/dx)  flows. 
Based  on  the  boundary  layer  characteristics  upstream  of  the  ramp  leading  edge, 
the  pressure  gradient  parameter  *k0  = ( *j"/Tw)  (dp/dx)  was  nominally  0.4  and  1.85 
for  Ramps  1 and  3,  respectively.  The  ramp  contours  and  a sketch  of  the  ramp 
installed  in  the  wind  tunnel  test  section  are  indicated  in  Figure  1.  Photo- 
graph of  the  ramps  installed  in  the  tunnel  is  shown  in  Figure  2. 

The  ramp  models  were  designed  to  replace  the  floor  plate  of  the  wind  tunnel 
test  section  and  to  provide  a smooch  continuation  of  the  lower  nozzle  block. 

In  order  to  avoid  flow  breakdown  in  the  test  section,  as  a consequence  of 
blockage  introduced  by  the  presence  of  the  models,  the  maximum  height  of  the 
ramps  was  restricted  to  1.8  cm  (see  Figure  1).  As  a result,  this  limited  the 
length  of  Ramp  3 to  about  9 cm,  although  this  was  found  to  be  adequate  for  the 
purposes  of  the  present  experiment. 

For  Ramp  1,  the  curvature  of  the  surface  was  started  at  the  leading  edge  of 
the  model  which,  in  turn,  mated  to  the  trailing  edge  of  the  nozzle  block. 
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\l  though  . t i 1 It  r material  v. . • > c.vd  to  smooth  tht  transition,  tlu  combined 
c 1 1 1*  v t.-  itirv.iturt  an!  unit  ; i-  I i. nncnt  resulted  it.  tht  generation  ol  a 

weak  slnn.  di  :•  t vr  nance  at  the  It  iclin  fd.'t  oi  tht  ramp.  This  is  apparent  in 
the  Spar.--- schl  lert  n photograph,  ol  Figure  la,  which,  otherwise  indicates  a 
disturbance- t ret  flow,  inct  st  itie  pr>  ssurt  measurements  indicated  that  tht 
shot k induced  pressure  risi  was  less  than  10  ol  the  total  pressure  increase 
along  the  ramp,  the  ellect  at  tht  leading  edge  shock  on  the  subsequent  de- 
velopment oi  the  i 1 ow  was  considered  negligible,  and  the  ramp  was  judged 
adequate  tor  use.  A different  strategy  was  adopted  in  tht  labrication  ol 
Ramp  3.  For  this  model,  surface  curvature  was  not  initiated  until  2.54  cm 
downstream  ol  the  leading  edge  with  the  portion  in  between  machined  llat. 

This  separated  the  curved  section  of  the  ramp  t rom  the  junction  with  the  nozitlt 
block  trailing  edge.  With  a tiller  material  again  used  to  smooth  the  transi- 
tion from  the  noggle  to  the  ramp  model,  it  was  possible  to  produce  a shock- 
free  flow  over  the  initial  portion  ot  the  ramp.  This  is  illustrated  in  the 
Schlieren  photograph  ol  Figure  3.'  , The  density  disturbances  which  are  apparent 
in  the  photograph  were  demons  era  tod  ov  pressure  measurements  to  represent  a 
continuous  compression  lan,  Furthermore , tlu  pressure  measurements  show  chat 
the  oblique  shock  observed  in  the  downstream  1 1 ow  originates  2.3-3  cm  down- 
stream ot  the  beginning  ot  the  curved  vamp.  The  shock  is  a consequence  oi 
coalescence  ol  pressure  waves  generated  t>v  the  ramp,  and  its  existence,  location 
and  origin  are  predicted  by  the  Method  ol  Characteristics  solution  to  the  tlow 
field. 

Each  ramp  was  provided  with  0.0t>4  cm  diameter  pressure  ports  aligned  along  the 
ramp  centerline  at  1.27  cm  intervals.  The  measured  surlace  pressure  dis- 
tributions tor  the  two  ramp  models  are  shown  in  Figure  4,  which  also  includes 
data  obtained  Iron  the  static  pressure  surveys  discussed  later.  The  observed 
scatter  in  t lie  data  can  he  attributed  io  small  inaccuracies  in  the  pressure 
measurements  and  to  slight  imperfections  in  the  surlace  contour.  In  both 
cases,  however,  the  measured  pressure  increase  is  sufficiently  linear  that 
the  f 1 ow  can  he  considered  as  characterised  by  a constant  pressure  gradient 
with  dp/dx  equal  to  1.2  mmllg/cm  and  5 mmHg/cm  lor  Ramps  1 and  3,  respectively. 
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PRELIMINARY  FLOW  FIELD  SURVEYS 


Prior  to  the  conduct  of  the  final  measurements,  qualitative  pressure  surveys 
were  carried  out  to  insure  two  dimensional,  di sturbance- f ree  flow  over  the 
ramps.  For  Ramp  1,  pitot  pressure  surveys  to  examine  the  flow  volume  over 
the  ramp  were  made  at  seven  axial  stations  ranging  trom  1.2  cm  ahead  of  the 
leading  edge  to  2.5  cm  upstream  of  the  trailing  edge.  At  each  station, 
continuous  pitot  pressure  surveys  extending  laterally  (z)  to  1 cm  either  side 
of  the  ramp  centerline  were  made  at  approximately  one  dozen  vertical  (v  * ) 
locations  from  the  surface  to  about  1.8  cm  above  the  surface.  In  addition, 
at  the  same  x stations,  continuous  pt  versus  y*  profiles  were  traced  at  five 
z locations  (approximately  = 1.0,  .4  and  0 cm).  Typical  records,  shown  in 

Figures  5 and  6,  indicate  the  flow  is  free  of  gross  disturbances  with  only 
marginal  cross  flow  effects.  A graphic  representation  of  the  flow  field  is 
presented  in  Figure  7 where  continuous  pt  versus  y'  profiles  obtained  at  .63  cm 
intervals  along  the  ramp  centerline  are  plotted.  Flow  is  from  right  to  left 
with  the  right  hand  profile  located  .63  cm  downstream  of  the  leading  edge  of 
the  ramp  and  the  left  hand  profile  corresponding  to  14.6  cm  downstream  of  the 
leading  edge.  The  leading  edge  shock  (actually  two  closely  spaced  weak  shock 
waves  also  visible  in  the  Schlieren  record  of  Figure  3a)  is  clearly  indicated. 

Similar  measurements  were  carried  out  for  Ramp  3 and  typical  results  are  shown 
in  Figures  8-10.  Figure  8 presents  lateral  surveys  of  static  pressure  at 
several  positions  above  the  ramp  surface,  while  lateral  traverses  of  pitot 
pressure  are  indicated  in  Figure  9.  Roth  figures  indicate  that  the  flow  is 
two-dimensional  and  relatively  disturbance- free  and  Figure  8,  in  particular, 
implies  no  cross  flow  near  the  center  of  the  ramp  model.  Figure  10  is  a 
pitot  pressure  map,  similar  to  Figure  7,  which  depicts  the  development  of  the 
flow  field  along  the  length  of  the  ramp.  Although  the  crossing  of  the  pitot 
pressure  traces  slightly  obscures  the  clarity  of  the  figure,  it  is  still 
possible  to  detect  the  oblique  shock  formed  by  the  coalescence  of  pressure 
waves  generated  by  the  curved  surface. 
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INSTRl'MENTATION 


Measurements  were  made  ot  mean  profiles  of  pitot  pressure,  static  pressure, 
and  recovery  temperature  across  the  boundary  layer.  Each  probe  was  mounted 
in  the  tunnel  separately  to  avoid  am  possibility  of  mutual  interference  be- 
tween probes  and  to  minimize  blockage  caused  by  the  presence  ot  the  probes. 

The  pitot  pressure  probe  consisted  ol  a .0152  cm  O.D.  cube  which  was  acid 
etched  to  provide  a gradual  taper  along  its  tip  to  its  .0076  cm  l.D.  The 
probe  was  attached  to  an  aerodynamic  strut  which  was  soldered  to  a conical 
body,  housing  a miniaturised  Kulitc  pressure  transduce!  (Model  VQH-250-10), 
located  at  the  front  of  a remotel\  driven  actuator.  Tht  c lost  coupling  be- 
tween the  transducer  and  the  probe  tip  afforded  . significant  improvement  in 
the  response  time  of  the  measurement . Tlu  transducer  sensitivity  of  approx- 
imately 10  tnmHg/mv  was  checked  by  t 'libration  prior  to  each  run.  The  estimated 
maximum  error  in  the  pitot  pressuit  measurements  we*  c t when  tin  probe  was 
located  adjacent  to  the  wall  nci  t!  * 1 v ov : f • ;u  i am}  and  diminished  to  0,25.. 

in  the  trecstrear  ovet  the  rear  ot  t hi  *.  s . V irious  probt  corrections  due  to 
rarefaction  and  viscous  « ttect  wen  found  • > i iu  li.i:  li  although  slight 
interference  iifeits,  cause  ! \ proxi:*it\  oi  <e,  wcii  observed  for 

positions  vi  rv  c lost  to  tin  w ill.  !h  m if  tec  I how  wet  , wen  verv  sm.ill  and 

restricted  to,  at  most,  a tew  positions  and,  flin  ion  , i .*  correitions  were 
made . 

Static  pressure  m»  asur*  lent.**  wi  r«  cult  using  n op  i vr-c  v l itulci  shaped  probe 
with  a 0.051  cm  diamt  ai  . Four  statii  pn  * till  port  loc.ited  equallv  spaced 
iround  the  probi  peripherv  win  1 Ol  md  i 0.6'  * from  hi  probe  tip.  The 
probe  connected  to  Ovnisco  Mode  I Ti  AFT  S5-2  pressur*  transducer  located 
outside  • In  win*!  tutvu  1.  Tht  i : tin,  ted  maximur  errol  in  the  st  itii  pressure 
measurement  varied  t row  2r  with  the  probe  positioned  at  tin  beginning  ot  the 
ramp  to  0.67.  at  the  rear  ol  K imp  t.  Viscous  interac  tion  lorrections  tor  the 
static  prot'i  were  found  to  c negli  n li  for  the  Mach  numher-Revnolds  number 
conditions  at  flu  presen:  tests.  However,  the  static  pressure  proliles  in- 
dicated a small  (10  j intei  terenei  i I tec  t near  tlu  surt.iee  ol  the  ramp  which 
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extended  out  to  about  y'  = .12  cm.  Since  the  static  pressure  varied  nearly 
linearly  through  the  remainder  of  the  boundary  layer,  it  was  possible  to 
correct  for  the  interference  effect  by  extrapolating  the  static  pressure  pro- 
file to  the  wall.  This  permitted  a redundant  determination  of  the  surface 
pressure  distribution.  These  results  are  shown  in  Figure  4 where  they  are 
seen  to  provide  excellent  agreement  with  the  direct  measurements  of  wall 
pressure . 

For  both  the  static  pressure  and  the  pitot  pressure  measurements,  particular 
care  was  exercised  to  insure  that  the  probe  axis  was  parallel  to  the  local  ramp 
surface  prior  to  the  survey.  Both  probes  were  found  to  be  insensitive  to  yaw 
for  yaw  angles  less  than  5 degrees.  Since  this  condition  was  easily  satisfied 
across  the  boundary  layer,  the  boundary  layer  pressure  measurements  were  un- 
affected by  yaw. 

A bare  wire  Ch-Al  thermocouple,  installed  in  a ceramic  tube  and  also  attached 
to  an  aerodynamic  strut,  was  used  for  the  total  temperature  probe.  The  thermo- 
couple wires  were  welded  to  form  a disc  0.013  cm  thick,  thereby  providing 
almost  the  same  resolution  as  the  pitot  pressure  probe.  The  probe  was  cali- 
brated in  the  wind  tunnel  freestream  to  determine  its  recovery  factor  versus 
Reynolds  number  characteristic,  which  can  be  expressed  as: 


T -T 
meas 

*=  0.9151  + .4799  x 10  (Re  n)  - . 02302  x 10  Re  „ 

o,u  o , D 

o 

The  probe  Reynolds  number  is  evaluated  at  the  total  temperature  Tq  using  a 
characteristic  length  equal  to  the  disc  thickness.  The  output  of  the  thermo- 
couple could  be  read  with  a resolution  of  0.01  mv,  corresponding  to  0.25rK. 

The  tunnel  stagnation  pressure  was  measured  with  a 0-S00  mmHg  Heise  pressure 
gauge  with  a least  count  of  1 mmHg  and  stagnation  temperature  was  sensed  by  a 
Precision  Digital  Temperature  Indicator  which  read  directly  in  degrees 


7 


d 


Fahrenheit  with  a resolution  of  Is  F (,b~K).  Although  the  ramps  were  not  in- 
strumented to  measure  surface  temperature,  on  the  basis  of  a similar  boundary 
layer  study  (3)  performed  in  the  same  facility,  it  was  assumed  that  T = 

.‘(45  To  = 30fh'  K which  corresponds  essentially  to  the  adiabatic  wall  condition. 
Measurement  of  surface  pressure  was  made  with  the  model  IT  APT  85-3  Dvnisco 
transducer  used  for  the  static  pressure  surveys. 

Measurement  ol  wall  shear  stress  was  made  using  Preston  tubes  with  0.1,  0.163, 
and  0.236  cm  O.D.s.  The  smallest  and  largest  tubes  were  sized  using  criteria 
available  in  the  literature  for  adverse  pressure  gradient  flows  lo  determine 
the  minimum  and  maximum  probe  diameters.  Again,  the  Dvnisco  pressure  trans- 
ducer used  tor  the  static  pressure  measurements  was  used  to  acquire  the 
Preston  tube  data.  Since  the  results  were  found  to  agree  within  5T.,  only  the 
data  lor  the  0.1  cm  diameter  probe  are  discussed  later. 

Photographs  of  the  several  probes  described  above  are  shown  in  Figures  11  and 
12.  In  each  case,  flow  is  from  left  to  right  and  the  probe  support  is  designed 
so  that  aerodynamic  disturbances  art  swept  downstream  of  the  probe  tip. 

5 . TEST  PROCEDURE 


Profile  data  was  acquired  by  first  locating  the  probe  adjacent  to  the  ramp 
surface  where  the  surface  was  located  using  a 10-power  microscope  with  a 
calibrated  grnticle.  For  ihe  pitot  pressure  profile,  the  probe  was  then  moved 
vertically  upward  at  selected  intervals  until  the  transducer  response  was 
sufficiently  rapid  to  move  the  probe  at  a constant  slow  rate.  For  the  other 
parameters,  the  variations  were  sufficiently  srnal1  to  permit  continuous  tra- 
versing across  the  entire  boundary  layer.  A voltage  signal  proportional  to 
probe  position  and  the  sensor  signal  were  ted  to  an  A/P  system  whose  output 
' >s  recorded  on  tape  cassette  to  form  a permanent  data  file.  Probe  position 
intervals  ranging  trom  0.5  mils  to  5 mi l (depending  on  the  sensor  and  the 
rate  of  change  of  the  measured  variable)  were  used.  The  recorded  data  was 
subsequent  1 v stored  in  the  Company's  main  computer  where  it  was  processed  via 
a time  share  terminal. 
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Although  the  probe  actuator  is  provided  with  two  degrees  of  freedom,  it  is 
constrained  to  move  in  a vertical  direction  (normal  to  the  tunnel  centerline) 
as  opposed  to  normal  to  the  surface.  While  the  probe  could  be  moved  in  both 
the  x and  y'  directions  in  order  to  track  the  normal  to  the  surface,  it  was 
considered  more  convenient  to  obtain  the  surveys  along  the  vertical  and  use 
a simple  computer  programmed  interpolation  routine  to  convert  the  data  to 
profiles  perpendicular  to  the  model  surface. 

6.  TEST  MATRIX 


All  tests  were  conducted  at  the  tunnel  stagnation  conditions  listed  in  Section 
2.  Mean  flow  surveys  were  conducted  at  the  axial  stations  shown  in  Figure  13, 
which  also  includes  the  location  of  the  wall  pressure  ports  and  illustrates  the 
start  of  the  ramp  surface  relative  to  the  trailing  edge  of  the  nozzle  block. 


SECTION  III 


DATA  REDUCTION 


The  d.ic.i  reduction  procedure  was  programmed  lor  the  Company's  Honeywell  66/40 
digital  computer  and  all  data  processing  was  carried  out  via  a time  share 
terminal.  A schematic  oi  the  data  reduction  routine  is  shown  in  Figure  14. 

The  routine  is  actually  comprised  of  a number  of  sub-routines,  each  designed 
to  complete  a specific  calculation.  The  output  of  each  sub-routine  is  stored 
in  a DATA  FILE  which  is  used  as  input  to  subsequent  sub-routines  and  which  can 
oe  accessed  via  a graphics  terminal  to  provide  a hard  copy  graphical  representa- 
tion of  the  file  contents. 

The  pitot  pressure,  static  pressure  and  recovery  temperature  profile  data  are 
recorded  separately  on  tape  cassette  during  the  boundary  layer  survey.  The 
tape  cassettes  are  then  fed  to  the  computer  creating  three  data  files  for  each 
survey  station.  These  files  serve  as  input  to  Program  BLSURV2  which  performs 
three  functions.  First,  the  data  is  converted  from  "as  read"  units  to  physical 
units.  Second,  since  the  y'  positions  lor  the  static  pressure  and  recovery 
temperature  profiles  differs  from  those  for  the  pitot  pressure  survey,  the  data 
lor  the  former  is  interpolated  to  provide  static  pressure  and  recovery  temper- 
ature data  at  the  same  v'  locations  as  the  pitot  pressure.  Finally,  mean  flow 
properties  were  calculated  by  means  of  standard  gasdynamic  equations  using  an 
iterative  procedure  to  account  lor  the  calibrated  recovery  temperature  character- 
istics of  the  T probe.  The  resulting  boundary  layer  profiles  were  stored  in 
DATA  File  ZMFLXXXX,  where  Z denotes  the  ramp  and  XXXX  represents  the  x station. 

A typical  printout  from  BLSl'RVZ  is  shown  in  Table  1.  Similarly,  to  demonstrate 
the  density  oi  the  data  points  and  the  quality  of  the  measurements,  plots  of 
F./P,,  E./P  and  I'/T  versus  v'  are  shown  in  Figures  15-17,  respectively. 

t o s o -Of  r 

Profiles  of  1’  and  To  lor  the  same  x station  are  presented  in  Figures 
me as  me as 

18  and  ll),  respectively. 

For  a given  ramp,  the  Files  ZMFLXXXX  serve  as  input  to  PROGRAM  NEWT LON  which 
uses  a simplt  lour  point  interpolation  scheme  to  convert  the  profiles  measured 
along  tlu  vertical  to  protiles  along  the  normal  to  the  surface.  At  a given  x 
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station  the  interpolation  is  carried  out  using  profile  data  along  the  vertical 
at  that  station  together  with  vertical  profile  data  at  the  nearest  upstream 
station.  The  interpolated  flow  field,  which  includes  profiles  along  the  normal 
at  all  x stations,  is  stored  in  a single  DATAFILE  NEWFLOKZ  where  Z denotes  the 
ramp.  The  content  of  NEWFLOW  is  shown,  for  example,  in  Figures  20,  21  and  22 
where,  respectively,  u,  P and  pg  have  been  plotted  against  the  distance  normal 
to  the  surface,  y.  Similar  plots  for  Ramp  3 are  shown  in  Figures  23-25.  In 
addition,  normal  profile  data  at  each  x station  is  stored  in  separate  DATA  FILES 
ZNFLXXXX.  The  latter  is  used  as  input  to  PROGRAM  VCOLES  which  correlates  the 
experimental  data  with  Coles  "Law  of  the  Wake."  This  calculation  provides  the 
boundary  layer  thickness  5,  the  wall  shear  stress  ' , and  the  wake  parameter 
Furthermore,  the  transformed  velocity  profile  is  stored  in  DATA  FILE 
ZVCOXXXX  in  u+,  y+  coordinates  and  in  DATA  FILE  ZVDEFXXXX  in  velocity  deficit 
coordinates.  A sample  printout  of  the  output  of  PROGRAM  VCOLES  is  shown  in 
Table  2.  Using  the  value  of  5 provided  by  VCOLES  and  either  DATA  FILE 
ZNFLXXXX  or  NEWFLOWZ  as  input,  PROGRAM  TBLJNDIM  calculates  non-dimensional 
profiles  of  y/b  versus  u/ue,  o/oe,  etc.,  as  well  as  the  integral  properties  of 
the  boundary  layer.  A sample  printout  of  the  results  of  PROGRAM  TBLJNDIM  is 
presented  in  Table  3. 

2 

As  shown  later,  the  streamwise  derivatives  of  the  flow  variables  u,  Ou,  eu“, 

^nd  p are  needed  to  extract  the  turbulent  transport  properties  from  the  mean 
flow  data.  To  accomplish  this  DATA  FILE  NEWFLOWZ  is  input  to  PROGRAM  RAMSHER. 

An  x survey  station  is  selected  as  a reference  location  and  the  associated  v 
values  are  denoted  as  y . The  remaining  profiles  are  then  interpolated  to  find 
u(s,y  ),  pu(s,yo),  etc.  The  streamwise  distance  s at  y0  is  easily  determined 
knowing  the  axial  location  x of  the  survey  station  and  the  local  surface 
curvature.  For  each  yQ  the  flow  properties  are  curve- fit  by  the  method  of 
least  squares  to  the  expression: 

F(s,yo)  - Fj  + F2s  + F3s2 
2 

where  F * u,  pu,  Pu  , or  p.  With  F^,  F^  and  F^  determined,  the  streamwise 
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derivatives  are  Biven  by! 
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Once  the  derivatives  are  determined  lor  . given  proiile,  PROGRAM  RAMSHER 

computes  "/"  versus  v and  creates  DATA  FILER  ZSHERXXX,  storm  - v , and 

the  streamwise  derivatives,  and  ZEDYXXXX,  storing  v,  ",  . and  nu/jy.  File 

ZEDYXXXX  is  input  tinally  to  FROG RAM  EDDY  which  calculates  the  mixing  length 

and  eddv  viscosity  as  a function  ol  y.  PROGRAM  RAMSHER  can  also  he  easily 

moditied  to  create  a data  tile  storing  either  u,  , -u.  , etc.,  vs  s lor 

data  data 

selected  v or  u , . , -'u„ etc.,  versus  s.  Accessing  these  files 

•o  curvet it  curvet  it 

via  the  graphics  terminal  provides  a plot  oi  either  the  experimental  or  curve 
fitted  flow  field  which  is  convenient  lor  visually  assessing  the  quality  ol 
the  data  or  the  adequacy  of  the  curve  fit. 

A listing  ol  the  programs  used  in  the  data  reduction  can  he  found  in  Appendix 
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SECTION  IV 

EXPERIMENTAL  RESULTS 


1.  CHARACTERIZATION  OF  THE  PRESSURE  GRADIENT 

In  his  treatise  on  the  incompressible  turbulent  boundary  layer,  Clauser" 
concluded  that  the  proper  parameter  to  use  for  characterizing  equilibrium 
profiles  in  a flow  with  pressure  gradient  is  B = (f*/TW)dp/dx.  For  compress- 
ible flows,  Alber  and  Coats'’  suggested  replacing  ~ with  = ( 5 "/ ~w)dp/dx  and 
in  their  study  at  M.x  = 4,  Lewis,  et  alb  found  indeed  that  using  Ej  provides  an 
improved  correlation  with  the  low  speed  data.  For  the  present  experiment,  the 
variation  of  ?(.  with  axial  position  x is  plotted  in  Figure  26.  It  is  shown 
later  that  "w  increases  with  x much  faster  than  and,  since  dp/dx  is  a 
constant,  then,  as  indicated  in  Figure  26,  3^  decreases  in  the  axial  direction, 
particularly  for  RAMP  3.  Figure  26  includes  the  data  of  Sturek  and  Danberg 
which  also  obtained  for  constant  dp/dx  using  a curved  ramp  similar  to  that 
involved  in  the  present  tests.  While  a direct  comparison  cannot  be  made  since 
the  axial  coordinate  x has  not  been  properly  normalized,  their  data  also  show 
that  Sfc  decreases  with  increasing  x confirming  the  present  findings.  Further- 
more, the  magnitude  of  8^  in  their  experiments  indicates  that  the  relative 
influence  of  their  pressure  gradient  should  be  much  larger  than  in  the  present 
case.  This  point  is  addressed  later  in  this  report.  The  data  of  Lewis,  et  al, 
is  also  shown  in  Figure  26  to  demonstrate  the  similarity  in  the  magnitude  of 
the  pressure  gradient  parameter  3^  with  the  present  results.  Since  their  test 
was  carried  out  for  an  increasing  pressure  gradient,  dp/dx,  their  values  of 
increase  in  the  streamwise  direction.  However,  for  equilibrium  boundary  layers, 
this  should  be  immaterial  and  only  the  local  value  of  ?,  is  significant.  Ref- 
erence  to  the  results  of  Lewis,  et  al,  will  be  made  in  the  next  section  when 
the  present  data  is  compared  with  the  low  speed  correlations. 

2.  CORRELATION  OF  THE  VELOCITY  PROFILES 

It  has  become  common  practice  to  compare  experimentally  measured  veloci tv- 
profiles  to  a well-defined  law  (e.g,.  Coles  composite  "wall-wake"  correlation)^ 
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that  describes  the  behavior  of  an  equilibrium  turbulent  boundary  layer.  Using 
an  appropriate  transformation  to  convert  the  compressible  data  to  an  equivalent 
incompressible  form,  this  comparison  permits  an  assessment  of  the  quality  of 
the  data,  represents  a means  for  evaluating  the  characteristic  boundary  layer 
parameters  (e.g.,  6,  cf,  etc.),  and,  in  the  non-constant  pressure  case,  assists 
in  isolating  the  effects  of  pressure  gradient  on  the  development  of  the  bound- 
ary layer  and  provides  a basis  for  comparison  with  other  experiments.  The 
comparison  of  the  experimental  data  to  the  classical  boundary  layer  profile 
involves  curve-fitting  the  transformed  data  to  the  classical  profile,  while 
iterating  the  values  of  the  unknown  parameters  until  the  rms  deviation  of  the 
curve-fit  is  minimized.  Details  of  the  curve- fitting  procedure  and  a discussion 
of  the  results  are  presented  below. 

Starting  with  the  conventional  mixing  length  expression: 

T = Tw  = ol 2 (du/dy)2  1 

and  combining  with  the  Prandtl  hypothesis: 

l = xy  2 

yields  the  following  relation: 


du+  _ ! p_  ^ du  _ l dy 

pw  ut  = * y 

Integration  of  the  above  expression  gives 

u"*”  = — tn  y+  + C 

H 

where : 

y+  s yu  /v 

J J T W 


3 


4 
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Equation  4 is  the  conventional  "Law-of-the-U'all"  which,  following  Coles,  has 
been  replaced  by  the  more  general  "Law-of- the-Wake " formulation  of  the  mean 
velocity  profile,  i.e.,  by: 

1 *T 

u+  = - £n  y+  + C + - W ( y / t ) 5 

In  Equation  5 ~ is  a parameter  representing  the  strength  of  the  wake  component 

of  the  boundary  layer,  W is  coles  tabulated  wake  function  which  can  be  approxi- 
? 

mated  by  2 sin“  ("^y/Zb)  and  the  constants  and  C are  given  their  incompressible 
values  0.41  and  5.0,  respectively. 

Equation  5 contains  three  unknowns:  6,  u-  and  ~.  Substituting  the  edge 
conditions  into  Equation  5 yields: 


u+  = - In  + C + - — 

e x v x 

w 


which  can  be  used  to  express  rr  in  terms  of  u_  and  6.  This  reduces  the  number 
of  unknowns  in  Equation  5 to  two,  whose  values  are  adjusted  until  the  data 
fit  the  equation  such  that  the  rms  error  is  a minimum.  Data  near  the  wall  and 
near  the  edge  of  the  boundary  layer  are  excluded  from  the  curve-fit  and  only 
data  for  which  y+  £ 50  and  y/6  s 0.9  are  used  to  determine  the  boundary  layer 
parameters . 

Although  the  curve-fit  procedure  is  restricted  to  the  v^all-wake  region,  it  is 
instructive  to  compare  the  experimental  data  to  the  "universal"  velocity  profile 
across  the  entire  boundary  layer.  In  the  sub-layer  region,  the  velocity  profile 
is  commonly  expressed  as: 


The  transition  between  the  sublayer  and  wall  regions  of  the  boundary  layer  has 

9 10 

been  examined  by  Spalding  (and  later  by  Kleinstein  using  a more  formal 
approach)  who  suggests  that  the  velocity  profile  in  this  zone  can  be 
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described  bv : 


+ + 


= u + exp(-.C')  texp(  .u+)  - I (u  ) 


where 


i (u  + ) = 1 + U+  + - ( U+)-  + ~ (.U+)  ' + 77T  (..uV 

D •_  H 


8 


Note  that  Equation  8 reduces  to  Equation  7 as  .;u  -*  0 and  that  Equation  4 
is  recovered  when  ..u+  ■ > 1.0.  Equations  8 and  5 are  used  to  represent 
the  universal  velocity  protile  across  the  boundary  layer. 

Typical  plots  ot  the  experimental  velocity  proiiles  in  translormed  coordinates 
are  shown  in  Figures  27  and  2S  lor  Ramp  1 and  Ramp  3,  respectively.  The  values 
ol  u+  have  been  determined  directly  Irom  Equation  4.3  using  the  measured  density 
proiiles  to  carry  out  the  integration.  This  avoids  reliance  on  analytical 
transformations  and  their  approximations.  The  Van  Driest  Transformation,  for 
example,  which  has  been  successfully  used  lor  flat  plate  boundary  layers, 
assume  constant  pressure  and  relates  the  density  to  the  velocity  via  the 
Crocco  relation.  Equation  8 has  been  represented  by  dashed  lines.  For 
clarity,  a plot  ol  Equation  5 has  been  omitted  since  the  differences  between 
the  experimental  values  of  u+  and  the  theoretical  values  cannot  be  resolved 
within  the  scale  used  in  the  figure.  In  the  region  y • 50,  the  data  tend  to 

lie  above  the  theoretical  curve  (a  feature  common  to  a considerable  body  of 
experimental  data)  with  the  discrepancy  reaching  a maximum  in  the  range 
10  < v+  < 20.  This  portion  of  the  velocity  profile  has  been  replotted  to  a 
larger  scale  in  Figure  29  where  lor  each  ramp  data  irom  the  dp/dx  = 0 survey 
station  and  the  most  farther  downstream  station  on  the  ramp  have  been  included, 
together  with  the  theoretical  proiiles  given  by  Equations  4,  7,  and  8. 

The  dp/dx  = 0 data  lor  Ramp  1 is  seen  to  be  in  excellent  agreement  with  Equation 
8 while  for  Ramp  3,  the  data  lies  slightly  above  Equation  8.  In  addition, 
it  is  observed  that  when  dp/dx  > 0,  then  lor  both  ramps  there  is  a small  but 
definite  increase  in  the  discrepancy  between  the  experimental  u and  Equation 
8 in  the  streamwise  direction.  The  reason  lor  this  is  not  immediately 


lb 


r 
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obvious,  but  it  reflects  in  the  mixing  length  calculations  in  Section  IV-6, 
and  is  discuss  ed  further  in  Appendix  B. 

For  both  ramps,  a plot  of  the  wake  function  W,  for  the  same  survey  stations 

shown  in  Figures  27-29,  is  presented  in  Figure  30,  where  it  is  compared  to 

2 

Coles  approximation  2 sin“  (r~y/2$).  Although  the  agreement  between  the  ex- 
perimental and  theoretical  wake  function  is  considered  reasonable,  there  is 

2 

a systematic  increase  in  the  difference  W-2  sin“  (-rry/25)  in  the  streamwise 
direction,  with  the  sign  of  this  difference  changing  from  - to  + as  the  outer 
edge  of  the  boundary  layer  is  approached.  Furthermore,  in  contrast  to  the 
sine  function  which  vanishes  as  y/5  — 0,  W remains  finite  as  the  wall  is 
approached. 

In  Reference  8,  Coles  discusses  the  effect  of  pitot  probe  errors  on  velocity 
measurements  near  the  wall.  These  errors,  which  arise  from  a variety  of 
sources  and  include  probe  interference  effects,  uncertainty  in  probe  position, 
and  the  influence  of  locally  high  turbulence  levels,  are  difficult  to  diagnose 
and  to  correct.  This,  in  fact,  is  the  main  reason  why  the  curve-fitting 
procedure  is  restricted  to  data  for  which  y+  > 50.  For  the  present  tests,  the 
discrepancies  between  the  data  and  the  universal  correlation  as  illustrated  in 
Figure  29  and  for  particularly  y/6  <0.2  in  Figure  30,  are  similar  in  trend  and 

g 

magnitude  to  those  associated  with  the  numerous  experiments  examined  by  Coles1 
and  are  not  considered  unusual.  Consequently,  the  curve-fitting  procedure  is 
assumed  to  provide  an  accurate  determination  of  the  parameters  6,  uT  and  ~. 

A summary  of  the  boundary  layer  parameters  obtained  from  the  curve-fitting 
process  is  presented  in  Table  2. 

4 

Clauser  defines  an  effective  displacement  thickness  of  the  turbulent  boundary 
layer  in  terms  of  the  transformed  velocities  as: 

u * -S  (u  - U ) d(y/S)  9a 

o e 
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,ind  tht  corresponding’  shape  1 actor  as: 


G = (u  - u )“  d(v/_)  9h 

o e 

where  tor  constant  pressure  layers  G and  .J  ' have  the  values  6.8  and  3.6, 
respectively.  Written  in  the  velocity  detect  iorm  used  in  Equation  9, 
Equation  S becomes: 


+ + 

u - u 

e 


:/'■  - 


-(W-2) 
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Typical  plots  of  u - u^  versus  y/_  for  survey  stations  located  just  upstream 
of,  at  the  midpoint,  and  at  the  rear  of  the  ramp,  are  shown  in  Figure  31  for 
Ramp  1 and  Ramp  3 in  order  to  illustrate  the  effect  of  the  pressure  gradient 
on  the  shape  of  the  velocity  profile.  For  Ramp  1,  the  pressure  gradient  is 
relatively  weak  and  is  nearly  constant  and  the  velocity  profiles  are  in- 
dependent of  streamwise  location.  The  pressure  gradient  for  Ramp  3 is  stronger 
and  initially  the  velocity  profile  is  distorted  (compare  x = 5. OS  cm  station 
to  x = 0 cm  station).  However,  since  decreases  with  x,  the  effect  of  the 
pressure  gradient  diminishes  and  the  shape  of  the  velocity  profile  at  the 
rear  of  the  ramp  is  almost  identical  to  that  where  dp/dx  = 0. 


A plot  of  the  boundary  layer  thickness  ? versus  axial  position  x is  shown  in 
Figure  32  which  includes  values  ot  obtained  from  the  curve-fit  to  the  Law- 
of-the-Wake  correlation  and  those  determined  directly  from  the  measured  vel- 
ocity profile.  The  latter  were  evaluated  ironi  visual  inspection  ol  the  profile 
by  selecting  the  v location  where  the  boundary  layer  data  merged  with  the  data 
in  the  external  stream  and  denoting  this  position  and  the  associated  velocity 
as  the  edge  conditions.  The  thickness  was  defined  then  as  the  y position 
where  u = 0.995  u(  . With  the  exception  of  the  forward  portion  ot  Ramp  3,  the 
esults  derived  from  the  correlation  are  in  excellent  agreement  with  those 
obtained  1 rom  ttu  profile  data.  The  values  ot  used  subsequently  in  this 
report  are  those  determined  I rom  the  curve-fit  of  the  velocity  profile. 
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Figure  33  presents  a plot  of  the  wake  parameter  " versus  axial  station  x.  For 
both  ramps,  ~ jumps  to  a peak  value  then  similar  to  9^,  decreases  as  the  rear 
of  the  ramp  is  approached.  However,  while  is  a maximum  at  the  leading  edge 
of  the  ramp,  the  variation  of  - with  x in  this  region  is  much  slower,  with  - 
not  reaching  its  maximum  value  until  x = 3 or  4 cm.  Since  the  pressure  waves 
generated  by  the  curved  ramp  surface  are  swept  downstream,  then  just  behind  the 
leading  edge  of  the  ramp  the  outer  portion  of  the  boundary  layer  retains  a 
memory  of  its  upstream  history  (i.e.,  the  flow  here  is  still  characteristic  of 
dp/dx  = 0).  This  is  also  the  wake  portion  of  the  boundary  layer  which  con- 
tributes largely  to  the  value  of  ~.  Therefore,  just  downstream  of  the  leading 
edge,  the  boundary  layer  is  not  in  equilibrium  with  the  local  value  of  Jj.. 
(Strictly  speaking,  the  procedure  for  curve-fitting  the  data  to  the  "Law-of- 
the-Wake"  is  not  valid  for  these  first  few  survey  stations  since  the  velocity 
correlation  is  restricted  to  equilibrium  flows).  This  is  indicated  in  Figure  34 
where  ~ is  plotted  versus  r^.  The  lack  of  correlation  between  ~ and  -|.  at  the 
first  two  ramp  survey  stations  for  Ramp  3 is  quite  apparent.  Figure  34  also 
includes  the  data  of  Sturek  & Danberg^  and  Lewis,  et  al^,  and  the  results  of  a 
number  of  low  speed  experiments  examined  by  Coles  & Hirst*'*'  and  correlated  by 
Lewis,  et  al^.  Interestingly  enough,  the  data  of  Sturek  and  Danberg^  lie  near 
the  upper  bound  to  the  spread  of  the  low  speed  data  while  those  of  the  present 
tests  fall  near  the  lower  bound  of  the  low  speed  data.  For  Ramp  1,  the  values 
of  5^  are  too  small  to  identify  any  specific  trend  of  r with  5u.  Finally,  a 

4 

plot  of  Clauser's  shape  factor  G versus  9^  is  shown  in  Figure  35.  Again,  the 
present  data  for  Ramp  3 lie  near  the  lower  bound  of  the  spread  of  the  low  speed 
data,  although  the  general  trend  of  increasing  G with  increasing  Sc  is  apparent. 
The  results  for  Ramp  1 lie  within  the  low  speed  data  spread  while,  for  compari- 
son, the  results  of  Lewis,  et  al  are  near  the  upper  bound  of  the  low  speed 

data  and  in  slightly  better  agreement  with  the  theoretical  results  of  Mellor 

12 

and  Gibson  , Excluding  the  flow  just  downstream  of  the  ramp  leading  edge, 
the  results  shown  in  Figures  34  and  35  are  in  agreement  with  the  earlier 
findings  of  Lewis,  et  al  ; namely,  that  the  boundary  layer  is  in  approximate 
local  equilibrium  throughout  the  adverse  pressure  gradient  region.  Thus, 
the  boundary  layer  profiles  are  characterized  by  local  conditions  only  and 
unaffected  by  the  fact  that  9^  is  not  constant. 
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A plot  ol  the  wall  shear  stress  r deterr, lined  iron  the  turve-lit  to  the  "Law- 

ol- the- Wake , " versus  axial  station  is  shown  in  Figure  36  which  includes  wall 

shear  measurements  made  with  the  0.1  cm  diameter  Preston  tuhe.  The  Preston 

13 

tune  data  was  reduced  using  the  Bra'dshaw-l'nsworth  correlation: 

r ")2 

— 35  9b  ■+•  oO  log,n(u  d/50'.  ) + 23.7  wlog,_(u  d/50'.  )_ 
w 

, 1 f n w 1 

+ 10"*>r  L (u  d/'  )L  ‘~b  - 2.0  J 

” W 

where  .p  is  the  Preston  tube  reading,  d is  the  tube  diameter  and  M = u_/a  . 
The  two  sets  ol  results  are  in  very  good  agreement  although  the  Preston  tube 
measurements  are  generally  lower  than  the  data  obtained  lrom  the  velocity 
correlation  with  a maximum  difference  oi  9\.  at  the  rear  ol  Kamp  3.  Both  re- 
sults, however,  show  a continuous  increase  in  - in  the  downstream  direction. 
To  be  consistent  with  the  selection  ol  the  boundary  layer  thickness  :,  the 
wall  shear  stress  determined  from,  the  curve-fit  to  Equation  5 is  used  sub- 
sequently in  this  report.  A plot  oi  the  local  skin  friction  coefficient  cj 
versus  Re-,  is  presented  in  Figure  37  which  includes,  for  comparison,  the  skin 
friction  coefficient  for  dp/dx  = 0 calculated  from  the  Karman-Schoenberg 

equation  together  with  the  Van  Driest  transformation  as  outlined  in  Hopkins 
14 

and  Inouye  . For  a given  Mc  and  Re.;,  the  maximum  deviation  between  the 
measured  c£  and  that  calculated  for  dp/dx  = 0 is  about  1CT..  with  the  measured 
values  somewhat  larger.  This  contrasts  with  earlier  findings  which  indicate 
that  Cf  decreases  as  is  increased.  Nevertheless,  it  appears  that  the  in- 
crease in  the  local  wall  shear  stress  is  the  consequence  ol  increases  in  the 
1 ocal  dynamic  pressure  brought  about  by  the  pressure  gradient, 

3.  NON-DIMENSIONAL  BOUNDARY  LAYER  PROFILES 

As  indicated  in  Section  3,  the  boundary  layer  profiles  measured  along  the 
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vertical  to  the  surface  of  the  ramp  (i.e.,  perpendicular  to  the  tunnel  axis) 
were  interpolated  to  provide  profiles  normal  to  the  surface.  With  the  boundary 
layer  thickness  f1  now  determined  as  described  in  Section  IV-2,  the  boundary 
layer  profiles  along  the  normal  can  be  non-dimensionalized.  At  each  survey 
station,  the  velocity  u at  the  position  y - 5 was  arbitrarily  assumed  equal  to 
0.995  ue,  providing  a means  for  specifying  ufi.  The  profile  data  was  searched 
to  find  y ■ yfi  at  u ■ ue,  and  the  edge  values  of  the  remaining  properties  were 
defined  as  their  values  at  ye>  The  edge  conditions  ug,  p and  Mg  and  the 
integral  properties  6*  and  0 and  listed  in  Table  2.  Because  of  the  variation 
of  static  pressure  across  the  boundary  layer  and  the  lack  of  a freestream  region 
of  uniform  flow,  it  was  necessary  to  modify  the  conventional  expressions  for  the 
integral  properties.  The  definitions  used  here  were  suggested  by  McLafferty  and 
Barber and  recommended  by  Sturek  and  Danberg^,  and  take  into  account  the  flux 
deficit  appearing  within  the  boundary  layer  referenced  to  "ideal"  properties 
calculated  with  the  experimental  static  pressure  profile.  The  integral  thick- 
nesses are  referenced  to  the  ideal  properties  at  the  wall  and  in  the  case  of 
constant  static  pressure  reduce  to  the  classical  definitions.  According  to 
this  interpretation,  the  integral  profiles  are  given  by: 

uw  6k  " Jo  (u  • u)dy 

PwV  6*  “ MS  ' pu)dy 

oo 

pw\s  " .ro  Pu(u,<  ■ u>dy 

where  the  ideal  properties  and  u^  are  calculated  using  the  measured  static 
pressure  profile  assuming  constant  total  temperature  equal  to  the  freestream 
T^  value  and  constant  stagnation  pressure  equal  to  the  test  section  total 
pressure  pQ. 

Profiles  of  u/ue  versus  y/6  at  each  survey  station  along  the  curved  surtace 
are  shown  in  Figures  38  and  39,  for  Ramps  1 and  3,  respectively.  Similarly, 


profiles  ot  Mach  number,  M,  versus  v/  tor  R.imps  1 and  3 are  shown,  respectively, 
in  Figures  -«0  and  41,  These  figures  indicate  that  for  each  ramp,  the  development 
of  the  boundary  layer  is  gradual  and  continuous  and  tree  ol  discontinuities.  For 
Ramp  3 in  particular,  tht  profile  plots  reveal  gradual  changes  in  the  sublayer 
thickness  which  occur  along  the  length  ot  the  ramp. 
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In  a recent  paper,  Whitfield  and  High  examined  the  etiect  oi  non-unitv  Prandtl 

number  on  the  total  temperature-velocity  relationship  lor  :ero  pressure  gradient 

flows.  The  classical  Crocco  relation  H = u/ut>  is  restricted  to  unity  Prandtl 

number,  an  assumption  made  to  eliminate  the  turbulent  shear  stress  terms  lrom 

the  combined  energy-momentum  equation.  As  a consequence,  the  Crocco  relations 

fail  to  predict  the  well-known  total  temperature  overshoot  observed  in  non- 

unitv  Prandtl  number,  adiabatic  boundary  layers.  To  overcome  this  problem, 
lb 

Whitfield  and  High  introduced  an  approximate  model  lor  the  turbulent  shear 
stress  distribution  and  derived  an  analytical  solution  which  provides  a reason- 
able .t^rt-emeni  with  experimental  observations  (e.g.,  see  Reference  3),  For  the 
present  experiments,  it  was  found  that  the  H versus  u/u(,  relation  is  insensi- 
tive to  the  pressure  gradient  although,  as  shown  later,  the  shear  stress  dis- 
tribution apparently  is  strongly  dependent  on  It  is  ol  interest,  therefore, 

to  examine  more  closely  the  Whitfield  and  High  1 solution.  The  combined  energy 
momentum  equation  can  he  expressed  as: 


d“h 

— r + 0 
du" 


1 dh 

P ) — +1'  (\-l)  M"  = 0 

rm  ,-  rm 

- du  di. 


II 


where,  using  the  Whi  1 1 i eld-High  notation,  tin  overbars  represent  normalisation 
with  respect  to  tht  freestream  values  (i.e.,  it  = h/hx,  u = u/u.o),  Whitfield 
and  High  assume  that  the  Reynolds  shear  stress  is  proportional  to  the  turbulent 
kinetic  energy  which,  on  the  basis  ol  an  earlier  studv^,  leads  to  the  lollow- 
ing  approximation: 


exp  ( - 4(y / ) 


-V 


- u 


) 


Using  che  definition  of  u , Equation  12  can  be  rewritten  as: 


— *=  — exp  ( - 4(y/5)5^“  ) 
w :w 
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With  the  further  assumption  that  the  velocity  profile  can  be  expressed  as  a 
power  law: 


u * (y/5) 


1/m 
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Equation  11  can  be  written  in  the  form: 


h — ^ - e (4  i u ' 1h+-y)-^+(l-  e)  A h = 0 
du“  d u d u 
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where : 


1 - Pr 


■Y  ■ 5/2  m 
A - (Y-l)  m£ 

Equation  15  is  a second  order,  non-linear  ordinary  differential  equation 
for  h(u)  which  was  solved  by  Whitfield  and  High  ^ assuming  a solution  of  the 
form: 


h(u)  = hQ(u)  + e hj (u)  + .... 

subject  to  the  boundary  conditions  for  the  adiabatic  case: 
dho(0)  dh^O) 

o m 0 

d u d u 


hou>  - 1 

h^l)  **  0 
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Returning  to  Equation  13  the  density  ratio  (using  a typical  measured 
adiabatic  wall  proiile  lor  >«■  = 3 and  dp/dx  = 0),  the  exponential  term  and 
the  ratio  “/ ~ have  oeen  plotted  in  Figure  42.  It  is  seen  that  the  density 
chance  across  the  sublayer  is  much  larger  than  the  change  in  exponential  tern 


so  that  ~/~v  reaches  a peak  value  of  almost  1.5  at  y/  ~ 0.15.  As  shown  later, 


the  shear  stress  distribution  in  Figure 


is  similar  to  that  for  it. 


0.* 


In 


lb 


a recent  study,  Sandborn  concluded  that  for  dp/d>;  = 0,  the  shear  stress 

distribution  ~/~  versus  y/i  is  essentially  independent  of  Mach  number  and  in- 
w 

sensitive  to  wall  temperature.  In  fact,  his  "oest  estimate"  of  versus 


y/~  is  closelv  approximated  bv  the  exponential  term  in  Equation  13  and 
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Figure  42.  Consequently,  the  Khitf ield-High  analysis  is  repeated  below  using 

the  following  expression  ior  the  turbulent  shear  stress: 


= u.  'v  exp  ( 


( - 


►<>7  05/:) 


Note  that  at  y/:  = 0 Equation  IS  is  equivalent  to  Equation  4.12, 
using  the  definition  of  u - , Equation  IS  becomes: 


IS 

Again , 


— = exp  ( - 4(y/5)J' ") 
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1 /m 

Assuming  u = (v/5)  ' , Equation  11  can  be  expressed  as: 


d“h  ^ j dh 

— — - f 4-u  — 

d“u 


+ (l-f)A  = 0 


20 


du 


where  Equation  20  is  now  a second  order,  linear  differential  equation  for 
h(u)  subject  to  the  boundary  conditions  given  by  Equation  17.  It  is  again 
assumed  that  the  solution  has  the  form  expressed  by  Equation  16  leading  to: 


A - -2-  A 

1 +§  -1-u  2 - € [f  -J.-2-  + 


4Ai 


• -T+2n  x 

1-u  J / 
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(->•+1)  (y+2) 


where  at  the  wall: 


. . A . A 

1 + o - c ( - 


4/Vy 


('Y+l)(ry+2) 


■ 


\ 
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Note  that  Equation  21  reduces  to  the  Crocco  relation  when  z = 0. 

In  the  present  nomenc lature , the  Crocco  parameter  H can  be  expressed  as; 

_ A _0 

h + ~ u“  - h 

2 w 

H = 

23 

1 + 4 - h 


so  that,  substituting  Equations  21  and  22  into  Equation  23,  we  have  for 
the  adiabatic  case: 

, 8 ^ u 


H = u* 


('.+1)(~, +2) 


(y+l)(y+2) 


Thus,  to  first  order,  the  H - u relationship  is  independent  of  both  Mcc  and 
Pr  although  it  does  depend  on  the  velocity  power  law  exponent  m.  It  is 
interesting  to  note  that  the  correction  introduced  by  retaining  the  shear 
stress  term  in  Equation  H represents  a departure  from  the  Walz  quadratic 
law  rather  than  the  linear  Crocco  relation. 


To  verify  that  the  power  series  solution  for  h (u)  converges  rapidly,  Equation 
20  was  solved  retaining  terms  of  the  order  of  e . While  the  resulting  solu- 
tion for  H now  shows  a dependence  on  it  differs  from  the  first  order 
solution  by  less  than  1%. 

For  the  case  of  constant  wall  temperature,  Equation  11  was  solved  subject  to 

the  boundary  conditions: 

h (0)  = h 

o w 

h j ( 0 ) * 0 


h (1)  = 1 
o 

hjd)  - 0 
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Following  the  same  procedure  used  tor  the  adiabatic  case,  we  obtain  lor  the 


u relation: 


H = u . 1 


(1  - G'+1  ) 


l+=7  - h 


l-+l)(v*-2) 


1 + ^ - h 1 - u 
2 w 


A / , - » V 

- ~ (1-u) 


Ir.  this  case  attain,  the  crocco  re  . tion  = u is  recovered  when  € = 0. 


The  consequences  of  the  prt  sent  solutions  tor  versus  u will  now  at  consid- 
ered. Fox  the  adi  t . . < is < . qu  tioi  . is  plotted  it  Figux  t • ; : ox  several 


values  ol  : and  is  compared  to  the  .."hitf  idd-Ki  .th  solution  (for  M* 


1 - . typi  l experiment  d (f or  th«  same  Mach  number) . It 


is  seer  that  shape  ol  ■ c b.  versus  u . urv<  _:veu  -\  Equation  2-.  is  quite 


sir.i  lat  to  ‘ t soluti.'-.  at  : tin  .s.  no 


(o.nd  to  the  data)  and  that  noth 


predict  tat  . 


ovei  • ' . how < : 


, Equation  2-1  predicts 


lam  r v Ives  v : . :.c  .i::nKi-:a  * solution.  In  addition,  increasing 


th(  \ • ol  if  I ion  - c losei  to  th«  da  t a . Ac  or-.par  i son  o : : t 

power  1 iw  vt  Lie  : : pr  otiit  t oi  several  values  ol  m to  a typical  experimental 


ter,'  pressurt  cradiert  , ;d; . tic  w 11  ve  loci  tv  proiile  is  shown  in  Figure 

This  plot  c:e:-onstr  • ::  t : iu  power  law  profile  does  not  provide  a good 

represent  ti  v ot  t ■ c;.  inasmuch  s the  value  ol  m which  tits  the  experi- 
mental protilt  inert,  st-  with  \ . In  Figure  A3,  this  implies  that  F.-u 


relation  stilts  tow  : t.  curves  wit:  increasing  values  ot  m as  u increases.  In 


view  ot  tin  set  sitivirv  o:  Equation  24  to  the  exponent  m it  would  be  ot 
. -tertst  to  solve  the  sic  t qt:  tior.  us  in.;  a more  realistic  ve  loci  tv  protile 


This,  howevtr.  would  pi o v r»  quir*  numei ical  solution 


It  should  t pointed  n.  • 1:.  tor  ; hi  di  tv, tit  cast  . the  temperature  ditlor- 


Tw  i-  not  1 u a ’in  tor  the  present  ease  whi  i a 


3 . i s on  1 y 


about  IS  - 20cC.  The  parameter  H is  quite  sensitive  and  for  example,  an  in- 
crement of  0.1  in  H represents  only  a 2°  change  in  the  local  total  temperature. 
Hence,  the  apparently  large  differences  in  Figure  43  correspond  to  only  a few 
degrees  in  absolute  temperature.  With  this  in  mind,  it  is  suggested  that  the 

T_  distribution  across  the  boundarv  Inver  is  insensitive  to  the  model  assumed 
0 

for  the  turbulent  shear  stresses.  This  may  explain  why  the  present  tests, 
where  the  pressure  gradient  produces  large  changes  in  the  shear  stress  dis- 
tribution, indicate  similar  results  for  the  variation  of  K versus  u.  This  is 
shown  in  Figure  45  when  H versus  u has  been  plotted  for  dp/dx  = C and  for 
forward  and  aft  position  on  both  ramps.  The  similarity  in  the  H profiles  is 
quite  obvious  and  in  general,  differences  from  the  dp/dx  = 0 case  cannot  be 
distinguished. 

Before  closing  this  discussion,  it  is  instructive  to  examine  the  solution  for 

the  constant  wall  temperature  case.  A plot  of  Equation  26  for  = 3, 

Frm  = ,8S,  m = 7 and  several  values  of  wall  temperature  is  presented  in 

Figure  46  where  it  is  compared  to  the  Crocco  relation,  the  NTiitf ield-High^° 

3 

solution  and  typical  experimental  data  . The  differences  between  Equation  26 
and  the  solution  of  Whitfield  and  High  are  quite  small,  both  are  in  good  agree- 
ment with  the  data,  and  it  is  apparent  that  even  for  modest  heat  transfer  rates 
(i.e.,  the  Tw/Toe  = .714  case)  the  classical  Crocco  relation  is  a good  ap- 
proximation to  tne  data. 

4 DETERMINATION  OF  STREAMW1SE  DERIVATIVES 

In  order  to  calculate  the  turbulent  shear  stress  distribution  at  survey  stations 
located  along  the  curved  ramp  surface,  it  is  necessary  to  first  calculate  the 
y profile  of  the  streamwise  gradient  of  ou,  bu“  and  p.  At  this  stage,  the 
flow  field  data  is  specified  in  terms  of  profiles  of  u,  t,  ru , etc.,  along  the 
normal  to  the  surface  at  several  streamwise  locations.  The  set  of  v locations 
generally  differ  at  each  survey  station.  Therefore,  a reference  station, 
generally  located  near  the  midpoint  of  the  ramp,  is  selected  and  the  v locations 
at  this  station  are  denoted  y . At  each  remaining  station,  the  profile  data  is 


27 


interpolated  to  determine  the  flow  properties  at  y . For  each  y then,  the 
properties  u,  tu , tu“  and  p are  curve-fit  to  a second  order  polynomial  in 
terms  of  the  streamwise  distance  s and  the  resulting  expression  can  oe  differ- 
entiated analytically  to  determine  d cu(yo, s)/ds , etc. 

An  illustration  of  the  newly  interpolated  flow  field  and  the  resulting  curve- 
fitted  flow  field  is  shown  for  Ramp  1 in  Figures  47-49  where,  respectively,  u, 
ou  and  p have  been  plotted  versus  s for  selected  y positions.  Similar  plots  of 
u,  cu  and  p versus  s for  Ramp  3 are  shown  in  Figures  50-52,  respectively.  In 
carrying  out  the  curve-fit,  particular  caution  must  be  exercised  concerning  the 
data  from  the  first  few  survey  stations  on  the  ramp.  The  pressure  waves 
generated  by  the  ramp  are  swept  downstream  and  the  leading  pressure  wave 
penetrates  the  outer  edge  of  the  boundary  layer  several  centimeters  downstream 
of  the  ramp  leading  edge.  In  this  region,  the  inner  portion  of  the  boundary 
layer  feels  the  influence  of  the  pressure  gradient  while  the  outer  portion  is 
characteristic  of  dp/dx  = 0.  Moving  downstream  from  the  leading  edge,  the 
inner  portion  influenced  by  dp/dx  > 0 grows  thicker  while  the  outer  portion, 
where  the  effect  of  dp/dx  = 0 persists,  tends  to  vanish.  Data  from  survey 
stations  within  this  region  should  be  excluded  from  the  curve-fit  since  the 
profiles  here  are  not  characteristic  of  the  equilibrated  boundary  layer  in  an 
adverse  pressure  gradient.  This  effect  is  illustrated  in  Figure  51  where  the 
data  at  x stations  0,  .63,  and  1.9,  particularly  at  the  outer  edge  of  the 
layer,  does  not  blend  smoothly  with  the  downstream  flow  field. 

A partial  test  of  the  validity  of  the  curve-fit  is  provided  by  examining  the 
changes  introduced  by  arbitrarily  eliminating  data  from  some  of  the  survey 
stations.  For  Ramp  1,  for  example,  the  results  obtained  by  retaining  the  rear 
survey  stations  and  eliminating  one  or  more  of  the  forward  survey  stations 
(starting  with  x = 2.74  cm)  differed  to  an  unacceptable  degree.  Only  a slight 
improvement  was  observed  by  reversing  the  procedure.  However,  a significant 
improvement  was  obtained  by  including  the  x = -1.27  station  in  the  curve-fit. 

In  this  case,  eliminating  downstream  stations  produces  only  small  changes  in 
the  curve-fit  parameters  and  provided  acceptable  shear  stress  distributions 
over  most  of  the  ramp.  Including  this  station,  although  it  is  located  upstream 
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of  Che  ramp,  was  justified  by  the  relatively  small  pressure  gradient  associated 
with  Ramp  1 and  the  fact  that  the  edge  properties  do  not  vary  significantly  in 
the  streamwise  direction  (see  Figures  47-49).  For  Ramp  3,  however,  it  was 
considered  necessary  to  exclude  the  data  from  the  first  three  ramp  stations 
from  the  curve-fit. 
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Illustrations  of  the  y variation  of  the  streamwise  derivatives  oou/cs,  6du“/cs, 
and  op/os  for  Ramp  3,  x = 7.62  cm  are  shown  in  Figures  54-56,  respectively. 
These  curves  are  not  considered  typical  since  the  magnitude  of  the  derivatives 
and  the  shape  of  the  curves  depend  on  s and  dp/dx,  but  they  demonstrate  that 
the  curve- fitting  procedure  yields  continuous  results  with  relatively  little 
scatter  in  the  derivatives. 


5.  TURBULENT  SHEAR  STRESS  DISTRIBUTION 
a.  Zero  Pressure  Gradient  Region 


Combining  the  continuity  and  streamwise  momentum  equations  for  a zero  pressure 
gradient,  adiabatic  boundary  layer  yields: 


o,2 

T—  T = , I Oil 

w ' o 3x 


iy  g 

(duz)  dy  - u .'o  (pu)  dy 


2 

Following  Sturek  we  assume  that  the  flow  is  locally  similar  (i.e.,  that  u/ue, 
o/oe>  etc. , are  functions  only  of  y/6)  so  that  Equation  4.27  can  be  rewritten 
as : 


w 1 d5 

2 

p u 5 dx 

CO  CO 


i y ou 
" 0 

P u 


u 

co 
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The  assumption  of  local  similarity  is  a convenient  approximation  and  is  not 
valid  near  the  wall.  However,  in  this  region  the  contribution  of  the  con- 
vective terms  is  quite  small  and  the  use  of  the  approximation  across  the 
entire  boundary  layer  is  justified.  With  this  assumption,  we  can  also  write: 

1 di  1 d > 1 cf 
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1 


— 


Substituting  this  relation  into  Equation  28,  wc  have: 


— “ 1 + x i,  f'  Pu^dy  - u ' pu  dy 
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which  can  be  evaluated  numerically  using  the  measured  profile  data. 


b.  Adverse  Pressure  Gradient  Region 


2 

According  to  Sturek*-,  the  equations  of  continuity  and  momentum  conservation 
for  a two  dimensional  boundary  layer  over  a surface  with  longitudinal  curvature 
are : 


Continuity 


|“  (pu)  + -^7  .(1  + ky) (pv  + p 1 v 1 ' 


)] 


1 


Momentum 


1 

1+ky 


du  , 
Pu  — + 
ds 


(Pv 


+ P ‘v  1 ) ^ + (pv  + 


P v 


u 


k 

l+k 
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1 ftp  . _ftT 
1+kv  fts  ftv 

Integrating  Equations  31  and  32  in  the  direction  y normal  to  the  surface 
and  combining  the  resulting  equations  yields  the  following  relation  for  the 
shear  stress  distribution: 


T " 1 + T 

W W 


i ^ ^ ->  I y ^ 

P (pu*”)  dv  - uP  , (pu)  dv 

O OS  * O 08 


(y  f fy 


. 2 


— (pu)  dv  uP  k dv  «f  p — dv 
o os  o ds 
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The  streamwise  derivatives  appearing  in  the  above  equation  were  determined 

using  the  curve-fit  procedure  described  in  the  previous  section.  The  values 

of  the  wall  stress  " were  taken  irom  the  correlation  of  the  measured  profile 
w r 

data  with  Coles  "Law-of- the-Wake . " 

c.  Results 


The  computed  normalised  shear  stress  distributions  for  Ramp  1 are  shown  in 
Figure  56  and  include  the  stress  distribution  in  the  zero  pressure  gradient 
region  1.27  cm  upstream  of  the  ramp  leading  edge  and  those  obtained  in  the 
adverse  pressure  gradient  region  at  stations  ranging  from  x = 5.27  to  12.79  cm. 
The  scatter  in  the  data  points  is  nil  and  for  clarity,  the  shear  stress  dis- 
tributions have  been  represented  by  continuous  lines  drawn  through  the  data 
points.  In  the  zero  pressure  gradient  case,  the  computed  shear  stress  dis- 
tribution is  in  good  agreement  with  expectations  based  on  Sandborn's  "best 
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estimate"  for  flat  plate  boundary  layers.  For  the  adverse  pressure  gradient 
region,  the  stress  distributions  indicate  a peak  at  y a 0.2  to  0.3  cm  (y/6  ^ 

.3  to  .45),  with  the  location  of  the  peak  value  shifting  away  from  the  wall  at 
the  downstream  locations.  In  addition,  the  shear  stress  remains  finite  at  the 
edge  of  the  layer  although  there  is  a systematic  shift  in  the  sign  of  the 
residual  " from  negative  to  slightly  positive  in  the  downstream  direction. 


If  the  data  input  to  Equation  33  were  completely  accurate,  the  shear  stress 
should  tend  to  zero  in  the  external  flow.  (Actually,  the  flow  in  the  external 
stream  is  not  uniform,  so  that  the  streamwise  derivatives  there  do  not 
identically  vanish  and  t remains  finite,  albeit  small).  Even  if  the  measured 
profiles  were  highly  accurate,  each  of  the  subsequent  manipulations  to  which 
the  data  is  subjected  introduces  an  uncertainty  which  reflects  in  the  final 
result.  Thus,  the  non-vanishing  stress  at  the  boundary  layer  edge  in  Figure  56 
is  believed  to  be  primarily  a consequence  of  the  data  processing.  It  is 
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possible  to  "correct"  the  data  so  that  the  shear  stress  does  vanish  at  the  edge 

*”) 

ot  the  boundary  layer.  For  example,  Sturek"  altered  his  prolile  ol  oou/os 
versus  y to  eliminate  a large  negative  residual  stress  at  y/'  = 1.0.  However, 
there  is  no  rationale  ior  changing  only  one  ol  the  derivative  terms  to  the 
exclusion  ol  the  others.  It  is  also  possible  to  alter  the  curve-lit  described 
in  Section  1V-4  by  excluding  some  ol  the  profile  stations.  Here  again,  there 
is  no  basis  lor  culling  the  data  used  in  the  curve-lit  except  lor  data  just 
behind  the  leading  edge  as  discussed  in  the  previous  section.  Consequently, 
the  data  has  been  left  unaltered  and  the  shear  stress  distributions  at  x 
stations  8.98,  10.25  and  11.52  cm,  where  */■  becomes  negligible  at  the  edge 
ol  the  boundary  layer,  are  considered  representative  ol  Ramp  1. 

Since  the  flow  properties  trom  y = 0 to  y = y all  contribute  to  the  magnitude 

of  " at  y = y,  any  correction  which  would  cause  " - 0 at  y = i would  also 

introduce  a change  ol  the  same  sign  and  a proportional  magnitude  to  the  peak 

shear  stress.  In  Figure  5b,  it  appears  that  il  ' was  adjusted  to  vanish  at 

the  edee  ol  the  lover,  then  the  stress  distributions  "/*  versus  v/  tend  to 

w 

approach  each  other.  This  implies  that  the  normalized  shear  stress  distribution 
is  insensitive  to  x location,  i.e.,  to  , and  is,  instead,  dependent  on  dp/dx. 
The  results  show  turther  that  even  a weak  pressure  gradient  * ^ 0.4  produces 

a peak  shear  stress  00“.  greater  than  the  wall  shear. 

The  shear  stress  distributions  lor  Ramp  3 are  shown  in  Figure  57  which  includes 
again  the  zero  pressure  gradient  result  and  the  results  lor  surveys  stations  in 
the  adverse  pressure  region  ranging  Iron  x = 3.18  to  7.b2  cm.  While  the  peak 
shear  stress  in  this  case  is  3 to  3.5  times  the  wall  value,  the  behavior  ol  the 
stress  distribution  as  a lunction  ol  x is  similar  to  that  observed  lor  Ramp  1 
and  most  of  the  comments  made  concerning  the  Ramp  1 results  apply  to  Ramp  3 
as  well.  In  particular,  the  shear  stress  does  not  completely  vanish  at  the 
edge  ol  the  boundary  layer  (although  t -*  0 in  the  external  stream).  The 
residual  shear  stress  at  y 5 is  generally  positive  with  a maximum  absolute 
value  of  0.4  "w,  which  is  probably  the  maximum  uncertainty  introduced  by  the 
data  processing  procedure,  and  is  now  only  107,  ol  the  maximum  shear  stress 
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in  the  boundary  layer.  This  data  shows  more  conclusively  that  the  normalized 
shear  stress  distribution  appears  to  be  dependent  on  dp/dx  rather  than  5^ 
which,  for  this  ramp,  decreases  by  a factor  of  almost  two  in  the  streamwise 
direction.  Sturek's^  results  indicate  a relaxation  effect  on  the  shear  stress 
profile  since  the  maximum  '■/"  at  his  forward  survey  station  (which  was  located 
near  the  mid-point  of  his  ramp)  was  one-half  the  peak  t/tw  observed  at  his 
downstream  stations.  However,  in  his  case,  dp/dx  varied  continuously  from  0 to 
a finite  value  at  the  forward  survey  station  where  it  remained  constant  over 
the  remainder  of  the  ramp.  This  contrasts  the  present  experiment  where  a 
constant  dp/dx  was  imposed  at  the  ramp  leading  edge  and  may  account  for  his 
observations.  It  should  also  be  noted  that  the  peak  _/"w  increased  quicklv 
downstream  of  his  first  survey  station. 

The  effect  of  pressure  gradient  on  the  turbulent  shear  stress  distribution 
across  the  boundary  layer  is  illustrated  in  Figure  58  where  stress  profiles 
t/t  versus  v for  dp/dx  = 0 and  representative  stations  of  Ramp  1 (dp/dx  = 

1.2  mmHg/cm)  and  Ramp  3 (dp/dx  = 5 mmHg/cm)  are  shown.  It  is  apparent  that 
the  stress  distribution  is  extremely  sensitive  to  dp/dx  with  the  peak  value 
of  ^/T  rising  from  60  to  300%  above  the  wall  value  for  the  weak  to  moderate 
pressure  gradients  used  in  the  present  tests.  It  should  be  emphasized  that 
while  t/tw  versus  y is  invariant  with  x location,  the  absolute  magnitude  of 
" is  increasing  since  rw  increases  with  x.  This  is  demonstrated  in  Figure  59 
where  t versus  y has  been  plotted  for  selected  x stations  along  Ramp  3.  Note 
that  the  maximum  shear  stress  r at  the  downstream  position  on  this  ramp  is  5-6 
times  greater  than  ^ in  the  dp/dx  = 0 flow  ahead  of  the  ramp. 

Since  3^  is  a decreasing  function  of  x for  a constant  pressure  gradient  flow, 
it  cannot  be  used  to  correlate  the  shear  stress  t/tw  which  appears  to  be  in- 
sensitive to  x.  However,  the  parameter  3^0,  where  and  are  evaluated  in 
the  zero  pressure  gradient  flow  just  upstream  of  the  curved  surface,  is  a 
characteristic  of  the  flow  field  and  remains  constant  for  constant  dp/dx.  A 
plot  of  the  peak  value  ol  t/tw  versus  8)<0,  including  data  from  Sturek's^ 
experiment  and  the  present  tests,  is  shown  in  Figure  60.  The  dashed  line 
representing  a linear  variation  of  ( j )m  wit'1  ^ko  wns  dr;,wn  through  the 
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most  representative  values  of  ('/''walPmax  obtained  from  cbe  present  tests. 
Although  the  data  do  not  follow  a linear  variation,  they  indicate  a consist- 
ent trend  and  demonstrate  a reasonable  agreement  between  the  results  of  the 
two  experiments. 


6.  TURBULENT  TRANSPORT  COEFFICIENTS 

With  the  turbulent  shear  stress  distribution  determined,  it  is  possible  to 
calculate  the  mixing  length  L,  and  the  eddy  viscosity  e,  using  the  following 
expressions : 


l 

~ 


T/P  u - 
e e 


1 h 


(D/Pe)  [ d(u/ue)/?(y/6)  ]2J 


34 


and 


(o/oe) 


, 2 

T/Peue 

d(u/ue)/d(y/6)  ] (5*/5) 


.35 


In  the  above  expressions,  t is  the  turbulent  shear  stress  and  is  obtained  by 
subtracting  from  the  stresses  computed  in  the  previous  section  the  laminar 
contribution  y,du/dy.  Because  of  the  very  high  density  of  data  points,  du/dy 
was  determined  directly  from  the  measured  velocity  profile  using  a simple 
differencing  scheme  and  no  attempt  was  made  to  smooth  either  the  velocity 
profile  or  the  variation  of  du/dy.  As  a result,  the  results  for  L/b  and 
e/ue6*  reflect  the  scatter  in  the  velocity  gradient  term.  A typical  plot  of 
d(u/ue)/d(y/5)  is  shown  in  Figure  61. 

Near  the  outer  edge  of  the  boundary  layer,  both  t and  du/dy  tend  toward  zero. 
As  a consequence,  large  errors  are  introduced  in  the  calculation  of  e/ue6* 
and,  particularly,  in  t/b.  For  this  reason,  when  y/6  > 0.8  these  quantities 
are  considered  unreliable  and  are  not  included  in  the  results. 
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A plot  of  ' /■  versus  v/5  lor  the  aero  pressure  gradient  boundary  layer  upstream 

oi  Ramp  3 is  shown  in  Figure  62.  The  scatter  in  the  data,  while  not  small, 

does  not  detract  from  a well-defined  trend.  Both  the  slope  ol  the  data  in  the 

wall  region  ( /i  = 0.4  y/5)  and  the  magnitude  of  '■  / - in  the  plateau  region  art 

in  good  agreement  with  the  conventionally  accepted  results  ol  Maise  and 

McDonald.-  The  influence  ol  the  adverse  pressure  gradient  on  ./'*  is  shown 

in  Figures  63  and  64  for  Ramp  3 and  Ramp  1,  respectively.  Because  of  the  scatter 

in  the  results,  the  variation  of  5/5  with  y / 5 has  been  represented  by  a curve 

drawn  through  the  mean  of  the  data  points  and  the  maximum  range  of  the  scatter  is 

denoted  bv  a vertical  bar  on  each  mean  curve.  Khile  ~/~  was  shown  to  he  in- 

w 

sensitive  to  x,  the  mixing  length  is  not  normalized  by  a wall  parameter  and, 
therefore,  -t/5  is  dependent  on  the  x station  and  increases  in  the  downstream 
direction.  From  Equations  34  and  35,  it  can  be  shown  that 

4/5  y <-/e)*  / (Su/&y) 


e 


u 

e 


a ( ~/c ) / (?u/jy) 


A close  examination  of  the  data  reveals  that  while  " and  0 both  increase  with 
x,  the  ratio  ~/o  also  increases  and  ou/ov  actually  decreases.  As  a consequence, 
for  a given  v/5  in  the  plateau  region,  both  l/i>  and  e/u  5’’^  increase  in  the 
downstream  as  shown  in  Figures  63  through  67.  Two  points  of  particular  inter- 
est are  apparent  in  Figures  63  and  64.  First,  in  the  adverse  pressure  gradient 
region,  the  slope  k (5/5  = k y/5)  in  the  wall  region  is  0.65  and  is  independent 
of  x and  dp/dx.  This  finding  is  identical  to  Sturek’s  observation,”  and  its 
implication  on  the  "wall-wake"  velocity  correlations  described  in  Section  1V-2 
is  pursued  further  in  Appendix  B.  Second,  the  magnitude  oi  f/5  is  similar  to 
the  values  found  by  Sturek,”  although  Sturek's  value  of  and  those  for 

Ramps  1 and  3 differ  by  as  much  as  a factor  of  9. 

A plot  of  the  normalized  eddy  viscosity  e/ugf’’^  versus  y/5  for  the  same  dp/dx  - 
0 station  shown  in  Figure  62  is  presented  in  Figure  65  where  it  is  found  to  be 
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in  excellent  agreement  with  the  universal 1\  accepted  results  ol  Maise  and 

20 

McDonald  ior  a similar  Mtt  and  Re..  The  eiiect  oi  adverse  pressure  gradient 
cm  the  eddy  viscosity  is  shown  in  Figures  ho  and  h7  where,  similar  to  the 
mixing  length,  the  data  points  have  been  replaced  by  continuous  curves  laired 
through  the  points  and  the  maximum  scatter  in  the  data  is  represented  by 
vertical  bars.  These  tigures  reilect  the  same  behavior  observed  lor  the  mix- 
ing length  and  show  that  the  eddy  viscosity  increases  with  x,with  somewhat 
larger  increases  apparent  ior  the  larger  value  ol  dp/dx.  In  fact,  the  maximum 

...  O 

value  ot  s/ue£ "k  is  even  larger  than  that  found  by  Sturek-  although  Sturek's 

value  ol  “ k0  is  twice  as  large  as  that  lor  Ramp  3.  However,  for  the  aero 

pressure  gradient  case,  Sturek  tound  e/u  5 to  be  only  hall  as  large  as  that 

•>0  e k 

predicted  by  Maise  and  McDonald-  and  this  discrepancy  may  also  be  reflected 
in  his  adverse  pressure  gradient  results. 
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SECTION  V 

SUMMARY  AND  CONCLUSIONS 


Measurements  have  been  made  of  mean  flow  profiles  at  several  streamwise 
locations  in  the  supersonic  turbulent  boundary  layer  over  a curved  ramp 
surface.  Two  ramp  models,  designed  to  produce  a constant  adverse  pressure 
gradient  flow,  were  used  with  7ko  s (dp/dx)(5^/"w)^  (with  (*£/'w)0  evaluated 
upstream  of  the  ramp  where  dp/dx  **  0)  equal  to  0.41  and  1.85.  Analysis  of 
the  profile  data  indicated  that: 

1)  With  an  appropriate  compressibility  transformation,  the  data  correlates 
with  the  well-defined  Coles  "wall-wake"  incompressible  velocity  profile. 

2)  Correlation  of  the  wake  parameter  - and  the  Clauser  shape  factor  C with 
the  local  pressure  gradient  parameter  rj.  is  in  agreement  with  the  low  speed 
data . 


3)  In  agreement  with  the  earlier  findings  of  Lewis,  et  al , ^ the  boundary 
layer  appears  to  be  in  a state  of  local  equilibrium  and  is  not  dependent  on 
upstream  history. 

4)  The  total  temperature  profile,  in  the  form  (TQ-Tw/Toe-Tw)versus  u/ue  is 
insensitive  to  the  pressure  gradient  and  is  similar  to  the  variation  for  a 
zero  pressure  gradient  boundary  layer. 

5)  The  skin  friction  coefficient  ct  was  lound  to  be  essentially  the  same  as 
for  dp/dx  = 0,  implying  that  the  observed  increases  in  wall  shear  are  a con- 
sequence of  the  increased  external  stream  dynamic  pressure  introduced  by  the 
pressure  gradient. 

Using  the  "indirect  method,"  the  flow  field  measurements  were  further  analyzed 
to  extract  the  turbulent  transport  terms  lrom  the  mean  flow  data.  Results 
show  that: 

1)  The  distribution  of  ~!~  , '1/6  and  e/u,,f,.  across  the  boundary  lover  for  the 

V v-  K 
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zero  pressure  gradient  data  of  the  present  experiment  are  in  good  agreement 

18  20 

with  the  earlier  findings  of  Sandborn  and  Maise  and  McDonald. 

2)  The  variation  of  the  turbulent  shear  stress  T with  distance  from  the 
surface  is  significantly  distorted  by  even  modest  values  of  dp/dx.  In  contrast 
to  the  zero  pressure  gradient  distribution,  when  dp/dx  > 0,  ^ increases  above 
its  wall  value,  reaching  a maximum  at  y/5  about  0.3  to  0.4. 

3)  The  normalized  shear  stress  distribution  "/tw  versus  y/5  is  independent 

of  the  local  3^,  although  the  peak  value  of  "/"w  appears  to  correlate  with  5,^. 

4)  The  maximum  values  of  -5/5  and  s/ug5^  for  the  adverse  pressure  gradient  flows 
reflect  the  increases  observed  in  the  maximum  values  of  ". 

5)  In  the  region  of  the  wall,  the  slope  constant  k in  the  expression  L/b  = 

k (y/5)  is  0.65  for  the  adverse  pressure  gradient  case  in  contrast  to  the 

zero  pressure  gradient  value  of  0.4.  The  value  of  k is  independent  of  5^  and 

2 

x and  is  identical  to  the  result  obtained  by  Sturek  for  5^Q  *=  3.45. 
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8.71  0 3 36 

1 0. 66667 

X. . 4 1 j ; I’i 

1 1 i'i  0 63 
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SUMMARY  OF  BOUNDARY  LAYER  PARAMETERS  - RAMP  NO 


RAMP  1 


Figure  2.  Photographs  of  Ramps  Installed  in  the  Wind  Tunnel 
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RAMP  1 

Figure  3.  Schlieren  Photograph  of  Ramp  Flow  Field 
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0.09  — O • ^WALL  ^EAS 

O ♦ PSTATIC  SURVEYS 


Fi^urt  4.  Surf -icc  Pressure  Distributions 


Z (cm) 


Figure  9.  Typical  Lateral  Picoc  Pressure  Surveys  at  Selected  v'  Locations 
Above  Surface  of  Ramp  3.  X Station  is  4.4  cm  from  Ramp  L.E. 
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A)  PITOT  PRESSURE  PROBE 


Photographs  of  the  Pitot  Pressure  and  Total  Temperature  Probes 
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typical 


Typical 


Figure  19.  Typical  Plot  of  TQ  versus  y1  (Ramp 

r me  as 


X STATION  (cm) 


Figure  20.  Profiles  of  u versus  Distance  y Normal  to  the  Surface  for  Ramp 


Figure  21.  Profiles  of  p versus  Distant 


X STATION  (cm) 


STATION  (cm) 


ante  y Normal  to  the  Surface  for  Ramp 


STUREK&  DANBERG  (7) 


gure  26.  Streamwise  Variation  of  Pressure  Gradient  Parameter 


Figure  30.  Distribution  of  Wake  Function  W Across  the  Boundary  Layer 


Figure  33.  Streamwise  Variation  of  the  "Wake  Strength  Parameter 


Figure  37.  Comparison  of  Skin  Friction  Coefficient  to  Zero  Pressure  Gradient  Results 
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Number  Profiles,  Ramp 


gure 


U/UE=  (Y/S)1^ 

DATA  FOR  dp/dx  = 0,  RAMP  1 


ssure  Gradient 


Figure  48.  Comparison  of  Experimental  Flowfield  to  Curve-fitted 
Flowfield,  Ramp  1.  Plot  of  Mass  Flux  versus  Curvili 
Distance  S at  Selected  y's 


Figure  j'/ . Comparison  of  Experimental  Flowfield  to  Curve-fitted  Flowfi 
Pnrnp  3.  Plot  of  Static.  Pressure  versus  Curvilinear  Distant1 


Figure  54 


tn.i  1 i zed  Turbulent  Shear  Stress  Distributions  Cor  Ramp 


Figure  57.  Normalized  Turbulent  Shear  Stress  Distributions  for  R imp 


of  Pressure  Gradient  on  Normalized  Turbule; 
Stress  Distribution  Across  the  Boundary  Lay< 


Figure  59 


Figure  60.  Variation  of  Peak  Shear  Stress  "/t  with  Pressure  Gradient 
Parameter  3^ 
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Eect  oi  Adverse  Pressure  Gradient  on  Eddy 
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APPENDIX  A 


DATA  REDUCTION  PROGRAMS 


This  section  contains  Program  Listings  for  the  various  computer  codes  used 
in  the  data  reduction.  Each  listing  contains  REMARK  statements  providing 
an  explanation  of  the  data  inputs  required  for  program  execution  and  in 
some  instances,  typical  input  is  illustrated  in  DATA  statements.  Figure  14 
of  the  text  indicates  the  sequence  in  which  the  programs  are  used  during  the 
overall  data  reduction  process. 
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PRECEDING  Page  hi  ANjK 


■malt* 


PROGRAM  BLSURV2 


10  PEN  BL— SUPVEY— £ 

11  PEH  P IS  THE  PRMP  NO.Q1  IS  730  MNH6.  Q£  IS  570  DE6  P.W  IS  THE  WRLL 

12  PEH  TEMP  540  DEG  P.S5  IS  THE  TO  PPOBE  DIR*. 005  IN.  HMD  N4  IS  THE 

13  PEH  NUMBER  DP  X PPOFILES.  X IS  THE  PPOFILE  STfKIN>  . P0»  R0»  BO  RPE  THE 

14  PEH  P0S<HHHG>  FOP  THE  PT.PS.TT  SURVEYS?  CO. DO. TO  RPE  THE  TO'S  FOP 

15  PEH  T0SUPVEYSJN1 .N2.N3  RPE  THE  HUMBER  OF  POINTS  IH  THE  PT.PS.TT  SURVEYS 

16  PEH  C1.C2  RPE  COEFF  IH  Y <CM> =C1*Y CMV)  +€£,  C3  IS  P <MM> /P <MV>  FOP 

17  PEH  PT  SURVEY?  C4.C5.C6  RPE  SRHE  FOP  PS  SUPVEY?  C7.C8  RPE  Y CONVERSION 

18  PEH  FOP  TT  SUPVEY  RND  C9.D1  RPE  FOP  T (DEG  F>  =C9VT  <MV>  *D1 . XI  <N>  RND 
18  PEH  S2<N>  RPE  X'IN)  RND  Pto/P0<X>. 

£0  FILES  F1?F£?F3?F4?F5 

£1  DIM  Y <£00>  .P'£00>  »Z<£00>  . S <£  0 0>  .W<£00>  »T<200> 

£5  DIN  XI  <I5>  » S2  <15> 

3 0 DIM  U C£  0 0>  . P 1 <2  0 0>  . 03  ( £ 0 0>  . T 1 <£ 0 0> 

35  PERD  X.PO.RO.BO.CO.DO.TO.N1. N£* N3 
40  PERD  P » G 1 . G£ » to. S5 > N4 
45  PERD  Cl .C2.C3.C4.C5.C6.C7.C8.C9.D1 
50  FOP  N=1  TO  N4 
55  PERD  XI <N>  » S2  <N> 

60  IF  X=X1  <N'>  GO  TO  70 
65  NEXT  N 

70  LET  S < 0>  =S£  <N>  *R 0 

75  LET  Z<0>=0 

76  LET  T < 0 > =W 

77  LET  to  < 0>  =0 

80  PRINT  " BDUNDRPY  LRYEP  SURVEYS.  X < CM> =“ X*£ . 54 

85  PRINT 
80  PRINT 
85  PPINT 

100  PPINT  TRB <15> . "PITOT  SURVEY" . "STRT I C SUPVEY" . "TEMP  SURVEY" 

105  PPINT 

110  PPINT  " P0<MMHG> ".PO.RO.BO 

115  PPINT  " T 0 ■■  DEG  K)  " . <C0-*-460)  ✓ 1 . 8.  <D0«-460> /\ . 8.  <T0*460> /\ . 8 

1£0  PPINT  " PW  • HHHG'>  “ » S£  <N>  ♦PO.  S£  <N>  ♦RO.  S2  CN>  *3*0 

1£5  PPINT 

130  PPINT 

135  PPINT 

140  PRINT  " I.IHRT  RPE  FILES  PTXX.YX.  PSXXXX.  TXXXX . PSTXXXX.  PLOPXXXX  “ 

145  INPUT  Fl*.F2*.F3f .F4S.F5S 
150  FILE  nl.Flf 
155  FILE  “2.F2S 
160  FILE  -3.F3* 

16J  FILE  “4.F4S 

170  FILE  -5. F5S 

175  PPINT 

180  PPINT 

185  PPINT 

180  PPINT 

195  PPINT 
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465 


PITOT  PRESSURE  PROFILE- 


PR  I NT  " 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT  " 

PRINT 

PRINT 

PRINT  "Y  <MV>  "*  “PT  ^NV>  “ , “YCCW>  “ » "PT  <'MWH6>  " • “PT/PO- 
PPINT 

FOR  L= 1 TO  N 1 
PERB  -1 » Y < L> > P < L ) 

PRINT  Y <L>  » P <L>  * 

LET  Y^L)=Cl*Y(L)t-C2 
LET  P <L> =C3*R <L> 

PRINT  Y <L>  * P <L>  * P <L>  /RO 
NEXT  L 


PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT  " 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT  “ 

PRINT 

PRINT 

PRINT  " Y <H W > “ PS (MV) “ - 
PRINT 

FOR  L=1  TO  N£ 

PERB  “£»  Z <!L)  » S <L) 

PRINT  Z <L> »S  <L>  » 

LET  Z <L> =C4*Z <L> +C5 
LET  S <L>  =C6*S  <L> 

PRINT  Z <L>  * S <L) » S <L> /RO 
NEXT  L 


STRTIC  PRESSURE  PROFILE- 


' Y <:CM>  "»  "PS(MMHG)  ' * "PS^PO-- 


PRINT 
PRINT 
PRINT 
PR  I NT 
PRINT 
PRINT  ” 

PRINT 
PRINT 
PRINT 
PRINT 
PRINT  ■ 

PRINT 

PRINT 

LET  T0=T 0+460 

PRINT  “Y  <f1V>  “ * "TT  <MV>  “ * "Y  <CN)  " » "TT  (DE6  K>  ” » “ HV T 0 - 


PECOVERY  TEMPERATURE  PROF  I LE " 


141 


r - r - O'  if.'  O'  o o c oj  oj  o~>  ri  rt  rt  \f‘  u"  x ■£■  r - r-  co  oj  j.  a.  o ® -«  oj  oj  cr«  oo  -<t  * in  y>  x ■•£>  r - r-  co  X ooc  c - — ■ oj  oj  o-.  co  n 
1 1 i 't  1 1 ii'  ui  y.  y>  y.  u"  ir>  u"'  y>  y.  y.  y*'  y-  iT'  y.  y.  y>  y.  m y.  -x  o x ■•£>  x x-  x x -x  x x x x i>  x x x x x 'Xt-i-r-f-t- 1-  i-  r - r. 


1 


PROFILE  DFtTB,  X=“X*€.54 


'Y  CCM>  ",  “PT-'PO"  , "PS  'PO",  "TTVTO' 


(Y  < J>  -Z  CK-l > > / CZ  <K ) -Z  CK  — 1 ) > 


FOP  L= 1 TO  N3 
PEftD  s3»WCL>  *TCL> 

PRINT  W<L>  * T <L>  » 

LET  kl  CL>  =C7*W  CL>  ♦C8 

LET  TCL>  = < <C9*T<L>«-I>l  >+46(» 

PRINT  WCL>  »TCL)/l.8,TCL>/T0 

NEXT  L 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PR  I NT  " ♦♦»»»»»»♦»♦♦»»»♦»♦•»♦»»»♦»»♦»»»♦■»♦»»»' 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT  “ CUNNRPY  OF  PROFILE  DFlTFt,  X=“X*£.54 

PRINT 

PRINT  “ Y <CH>  ",  "PT-'PO"  , "PS-'PO",  "TT/TO" 

PRINT 

SCPPTCH  s4 
FOP  J=1  TO  N 1 
FOP  k = l TO  N£ 

IF  Z CIO  > Y CJ>  GO  TO  60  0 
NEXT  K 

LET  03  CJ>  =S  <K— 1 > *<S  CIO -S  <fc-l>>  ♦<Y  < J)  -Z  <K-1)  ) ✓ CZ  <K>  -Z  <K-1>  > 

FOP  L=1  TD  N3 

IF  W CL>  > Y C J)  GO  TO  6£0 

NEXT  L 

LET  T 1 < J)  = T < L-l  > CT  CL>  — T CL— 1 > > ♦ CY  < J)  -W  CL-1 > > ' <W < L> -W CL- 1 > > 

PRINT  YCJ>  ,PCJ>/PO,Q3CJ)^ftO,Tl  CJ>/T0 

WRITE  ~4.YCJ>  ,PCJ>/PO,03CJ>/PO,T1  CJ>/T0 

NEXT  J 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

PRINT  " 

PPINT 

PRINT 

PPINT 

PPINT 

PRINT 

PPINT  " NEPN  FLOW  PROFILES  CTBL-J) » PRNP  NO"P 

PRINT 

PPINT 

PPINT 

PPINT 

PPINT 

PPINT  " P0,N'N2  ="01*133.3 

PPINT  " TO. DEG  ► ="Q£/1.8 

PPINT  " Tl.l.DEG  K ="W'1.8 

PRINT  ‘ P<:.MHW6  ="SP(  N>  1 


NEPN  FLOW  PROFILES  ■ TBL-J) » PRNP 


P0»«/H2 
TO- DEG  ► 
Tl.l.DEG  K 
P<C.  MMRR 


="01»133. 3 
= "02-'  1.8 
="W' 1 .8 
= " S?  ( N>  1 


745 

750 


76  0 
765 

77  0 
— 

78n 
785 
790 
7Q5 
8 0 0 
805 
81  0 
815 


PPIMT  " X.CM  = '»2.54 

PPIMT 

PPIMT 

PPIMT 

PPIMT 

PPIMT  "Y  ‘CM  - "M"  . “LUM/SEC"  . "PHD  ‘K6^M3'>  ” * “PE  <H-1> " 

PPIMT  TAP ■ 15> • “T  <TEG  K>  “ » “TO <DE6  K) “ * “PS  <MMHG> " * “PHDU  <KG^H£SEC) “ 
PPIMT 

SCPATCH  »5 
FDP  M= 1 TO  Ml 

LET  P2=  - P ‘M~-  *+40)  / < Q3  ♦A  O') 

IF  P2' 1.89886  GD  TD  880 

IF  P2  <1  GD  TD  965 

LET  Ml  =SOP  <"5*  ‘ < <P2)a(1/3.5>>-1>> 

GD  TD  860 


880  LET  Ml  = > .5*P2>'‘‘1/1.6> 

825  LET  Fl=»  1.8*m*fin~3.5 

830  LET  F2='6'  ■ ,7*M1*M1>-1>>"'£.5 

835  LET  P 1 =F 1 ♦F? 

840  LET  21=  - P2-P1> -'P£ 

845  IF  ABS  »"Z1>  <-.  001  GD  TD  860 

850  LET  M 1 =M 1 +Z 1 

855  GD  TD  825 

860  LET  T3=T  1 ‘M'- ♦O2/T0 

865  LET  F=l  + < . 2*+11^MO 

87  0 LET  Cl  = . 001623*G3‘N'<^F«01'' <T3^A0> 

875  LET  82=4*.  0 1 *411  ♦SQP  < T3/F> 

880  LET  73=.  0 000  0 00£27^<T3~1 .5> / <T3+199> 

885  LET  S4=S  1 ♦t 3 ♦35  - ‘ 1 £^S3> 

880  LET  76=.  915094+.  0004799*SOP  <24> 0000230237*S4 
885  LET  36=<S6*'  1 - ‘ 1 •- F)  > > ‘ 1 /F> 

800  LET  T 4= ■ T 1 1 M - ♦OE> / 1 T 0*S6) 

805  IF  HUS  - T3-T4‘ • =. 5 GD  TD  9£0 
810  LET  T 3=T 4 
81 5 GD  TD  865 

820  LET  PI ‘ M - =$  1 »3£ . £♦ 16. 04 
825  LET  T5=T4/F 

830  LET  Ll  ‘ M'*  =M  1 ^4 8 . 01  ♦SOP  ‘ T5 ^ ♦ . 3 048 
835  LET  Vl  = . 0000000227^‘ T5~l  . 5> • T5-M99> 

840  LET  77=  ‘ ‘PI  1 M ' ♦Li ‘M>  > / < 1 2^V  1^1  57 . 4£6>  > »1  00 -'2  • 54 
845  PPIMT  Y-N-'  .Ml  •U‘M'>  .PI  ‘M'‘  ,37 

850  PPIMT  TAB  - IS'1  . T5- 1 . 8 • T4/ 1 . 8 » G3 ‘ M ‘ ♦©  1 / A 0 • U - M''  ♦P 1 ‘ M1- 
855  PPIMT 

860  WPITE  ~5 • Y ‘ M > • U ‘ M>  • P 1 ‘ M - • U ‘ M ' •f*  1 < N>  * T4/ 1 . 8 
865  ME XT  M 
870  PPIMT 
975  PPIMT 
980  7TDP 
1000  DATA 
1010  DATA 
1020  DATA 
1030  DATA 
1040  DATA 
1050  DATA 
1060  DATA 

Aqqci  enp 


3.25*730*730.5*730.5. 1 10* 1 09, 1 10. 146. 132. 123 
3*730.570.540. . 005. 15 

. 00 020718*  1 . 026531  . . 09687.  . 00020718.  1 . 0944  <»*  .04113*.  00020718 
1 . 027314. . 44218.32.782 

— 1 . 25 * . 02986 * - . 75 * . 03071 » — . 25 • . 02975 » 0 * . 02986 • .25*  . 03212 
. 3*5-  * . 04  * 1 . 25  • • 04732  * 1 . 5- * . 0536 1*  1. 75  * . 05955  • 2 . . 065  04 
2.25* . 06899*  2 . 5 * . 07207 *2.75* . 07882*3* . 08277 *8.25* . 092 06 
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PROGRAM  NEWFLOW 


10  PEM  NEWFLOW 

11  PEM  01  I S THE  NUMBEP  OF  THE  FIPST  X STATION 

13  PEN  03  IS  THE  NUMBEP  OF  THE  LFlST  X STATION  MINUS  ONE 

13  PEM  B IS  THE  PAMP  NUMBEP 

14  PEM  X < N > IS  THE  LOCATION  OF  THE  X STATION  -INCHES') 

15  PEM  D-N'-  IS  THE  HEIGHT  OF  THE  PAMP  <CM>  AT  X 

16  PEM  A-N'>  IS  THE  SLOPE  OF  THE  PAMP  •PAD-'  AT  X 

17  PEM  ENTEP  01.G2.B  AT  LINE  1000 

18  PEM  ENTEP  1ST  5 INPUT  FILES  AT  303 

19  PEM  ENTEP  PEMAINING  INPUT  FILES  AT  304 

30  PEM  A TOTAL  OF  10  FILE  NAMES. PEAL  DP  FICT IC I DUS. MUST  BE  INPUT 
100  FILES  Fi;F3.*F3 

110  DIM  VI <1 18) .U1 <1 18) .PI <1 18) .PI <1 18) .T1 <118) 

130  DIM  Y2<1 18)  .U2<  1 18)  »P£U  18)  »P2<1 18)  »T2.1 18) 

135  DIM  A < 1 5)  » P 1 5)  » X < 1 5) 

130  PEAD  G1.03.B 
132  FOP  N=1  TO  0341 
140  PEAD  X <N) . D (N) » A <N) 

143  NEXT  N 

150  PPINT  " INTEPPOLATED  PAMP  FLOG  FIELD" 

160  PPINT  “ PAMP  NO"B 

170  PPINT 
180  PPINT 
190  PPINT 

300  PPINT  “ UHAT  APE  PPOFILE  INPUT  FILES" 

303  INPUT  Flf.F3S.F3f.F4*.F5S 

304  INPUT  F6S.F7S.F8S.F9f.F0S 

313  PPINT 

314  PPINT 
316  PPINT 

318  SCPATCH  -3 
33  0 FOP  0=171  TO  03 

333  ON  0 GD  TO  £70. 224 . 23 0. 236. 243 . 348 . 354 . 260. 366 

334  LET  F1S=F3S 
336  LET  F2S=F3S 
338  GO  TO  370 
330  LET  F 1 S=F3S 
332  LET  F2S=F4S 
334  GO  TO  370 
336  LET  F 1 S=F4S 
338  LET  F2S=F5S 
340  GO  TD  370 

3 '2  LET  F 1 S=F5S 
£44  LET  F3S=F6S 
346  GO  TO  370 
248  LET  F1S=F6S 
350  LET  F2S=F7S 

253  GO  TO  270 

254  LET  F 1 S=F7S 
256  LET  F2S=F8S 
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358  GO  TO  370 
360  LET  F1S=F8S 

363  LET  F3S=F«S 

364  GO  TO  370 
366  LET  F1$=F«S 
368  LET  F£*=Fn* 

370  FILE  “1.F1S 
380  FILE  -3.F3J 

373  PPINT  " X ■ CfT1  = “ X 1 0+ 1 ^ *3 .54 

374  PPINT 
376  PPINT 

378  PPINT  "V  'CPO  " . "U<M/SEO  " » ~PH0<KG'H3>  * , "P  " . "TO  - DEG  F>” 

383  PP I NT 

384  FOP  N= 1 TO  118 

390  PERB  si  • Y 1 f N)  • U 1 <'N>«P1  <N>  * P 1 * T 1 <N'1 

400  PERD  s8,Y3<M>  jUSW  »R3<N>  .P3<N>  . T3  <"«> 

410  NEXT  N 

413  LET  Y1  ( 0>  =0 

414  LET  Y3«'0>=0 
430  FOP  N= 1 TO  100 

430  LET  Y=Y3  <N>  *003  <R  < Q-M  > > 

435  LET  S='X'.G+l>-X<G>>»e.54“Y»SIN<RCQ-M>) 

437  LET  S 1 = <X  i O-M  > —X  <Q>  > ♦3.54 
440  LET  Z=Y*CO£  ■ R <Q«-1 ) > ♦D  (Q+l  > 

470  LET  N=N 

480  FOP  L= 1 TO  100 

490  IF  Y1 <L>  *CO£ (R  >0>  > =>Y  GO  TD  510 
500  NEXT  L 

510  LET  Zl=Y2<m-*-D<Q-*-l> 

530  LET  Z2=Y3<'N-n-*-B<Q+l'> 

530  LET  Z3=Y1 <L)4D<Q) 

540  LET  Z4=Y1 <L-l>^r<G) 

550  LET  Z5=Z3+  * Z 1 -Z3> ♦£/£ 1 

553  LET  Z6=Z4>*Z3-Z4')»Sx£l 

554  IF  Z6> Z GO  TO  56  0 
556  IF  Z5<Z  GO  TO  564 
558  GO  TD  570 

560  LET  L=L— 1 

563  GO  TO  530 

564  LET  L=L+ 1 
566  GO  TO  530 

570  LET  VI  =ll  1 <L— 1 > ■♦•  <U£  <M— 1 > — U1  <L  - 1)  ) ♦S/Sl 
58  0 LET  V3=U1  <L)  ♦ <U2  — Ll  1 <L>  ) ♦S/'Sl 

590  LET  U=V  1 ♦ < V3-V 1 ) ♦<  Z-Z6>  s <Z5— Z6> 

600  LET  V1=P1  <L-1>«-<P3  <H-1>-P1  <L-1>>»S/S1 
610  LET  V3=P  1 (L>*-<P2<'f1'>  — P 1 ^L)  > *S/S  1 
630  LET  P=V  1 ♦ CVS-V 1 ) ♦•<  Z-Z6)  ✓ < Z5— Z6> 

630  LET  V1=P1  <L-1>«-<P3  <H-1>-P1  (L-1>>*S/S1 
640  LET  V2=P1  (L)  <P2  <W> —P  I (L>>*S'S1 

650  LET  P=V  1 < V3-V 1 ) ♦< Z— Z6>  ✓ <Z5— Z6> 

660  LET  V1=T1 <L— 1 > ♦ < T2 (M- 1 > -T 1 <L- 1 > > ♦S^S 1 
670  LET  V£=T1  (L)+-<T2^H)-T1 
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680  LET  T=V1+<V8-V1>*<Z-Z6>' <Z5-Z6> 

690  PRINT  Y*U,R,P,T 

700  WRITE  «3, Y»U»R,P,T 

705  NEXT  N 

710  PRINT 

720  PRINT 

730  PRINT 

740  NEXT  0 

750  STOP 

1 000  DRTR  1,1,1 

1010  DRTR  1 • 077, . 07601 ,. 044472, 1 .577, . 14171 , 
1020  DRTR  2. 077, . 22501 0721 1 ,2. 535, . 31587, . 
1030  DRTR  3.  035*  . 43013, . 095679,  3. 535*  . 5“i94T, 
1040  DRTR  4. 035* .70324, . 1 1828,4.535, .861 18* . 
1050  DRTR  5. 035*1. 03308, . 1399, 5.535, 1 . 2 188Q» 
9999  END 
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. 058744 
083673 
.10717 
12915 
. 15066 


PROGRAM  VCOLES 


T 


10  PEW  VCOLES 

11  PEW  FIS  I S WFLOXXXi  FPOW  WEWFLOW> . F2S  IS  VPEXXX*  F3S  IS  VCDXXXX.  X IS 

12  PEW  X STR-IW''.  P IS  PW  ' WWH6)  * T IS  WALL  TEWP=540  PEG  P,  P IS  FIPST 

13  PEW  GUESS  FOP  PELTR<CH>»  C7  IS  FIRST  GUESS  FOP  TRUW'PSF'i.  W1  IS  1ST 

14  PEW  GUESS  FOP  WO  OF  Y POIWT  WHERE  Y*  =50.  W2  IS  FIPST  GUESS  FOP 

15  PEW  WO  OF  Y POIWT  MHEPE  Y'PELTR*. 9*  W3  IS  WO  OF  POIWTS  IW  WFLOXXXX. 

16  PEW  SI  RWP  B 1 RPE  cowvepgewce  cpitepir  fhp  truw  hwp  peltr. 

100  FILES  Fi;F£!F3 

110  PIW  LU200.)  .Y<200>  >Z(1  000)  »P1  <200>  .01  <200>  »T0<200> 

115  PIW  G'  1 00~>  . P < 1 00> 

116  PIW  U3C200) 

117  PIW  U5  < I 00' . U6  < 1 00> 

120  PERP  F1S.F2$.F3*.X.P.T.F.C7,W1.W2»N3 

122  PERP  SI. FI 

123  LET  ► 8=W2-M 

124  LET  K 9=W 1 

130  PPIWT  " COPPELRT I DW  OF  VELOCITY  PROFILE.  X STR  =“X»2.54 

132  PPIWT  " COLES  V^P  GEWEPRL I ZEP  COPPELRTIOW 

140  PPIWT 
150  PPIWT 

152  LET  F2  =. 00162*P-T 
165  LET  P7=0 
170  FILE  “ 1 » F 1 S 
175  FILE  S2.F2S 
ISO  FILE  -3.F3S 

181  LET  U'0>=0 

1 82  FOP  W= 1 TO  W3 

184  PERP  si . Y < W>  . LI  >:  W"1  « PI  < W"*  .01  <W">  • TO  <N> 

186  LET  Y >W' =Y < W> '2.54 

188  LET  Ll 1 W'1  =U  ' W>  ♦S . 28  08 

18^  LET  PI  «W>=P1  <W>/'516.488 

1Q0  LET  08=  1 SOP  <P  1 < W'1  'P2.'>  > ♦«  U 1 W'  -Li  <W- 1)  > 

182  LET  P7=P7^08 

183  LET  U3  1 W"*  =P7 
1 Q5  WEXT  W 

200  PPIWT 
210  PPIWT 

250  LET  v=.  0000000227»*  T~  1 . 5)  ✓ ■ • T-M*^  *P2'> 

281  LET  H=1 

282  LET  P'H^P-2.54 

285  LET  V5=  0 

286  LET  V6=  0 
280  LET  W= 1 
300  LET  S=C7 

310  LET  F = S0P'S^P2'> 

320  GOSUB  700 
350  LET  H=W-M 
360  LET  S=S«-S1»S 
370  LET  F = SPP  < S - PS'1 
380  GOSUF  700 

400  IF  Z<W'''Z'W-0  GO  TO  420 
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410  GO  TO  3?  0 
42  0 LET  M=M-*-l 
430  LET  S=S-S1*S 
440  LET  F=SQPiS/P2> 

450  GOSUB  700 

470  IF  Z<M>>Z<M-1)  GO  TO  500 
48  0 GO  TO  420 

490  PPIMT  " Y"  » " Y/D"  • "Y4-“  . “U-*-CRLC"  » “U+WEftS" 

491  PPIMT  TRB < 15> • 'Wr 
495  PPIMT 

500  LET  S=S«-S1*S 
505  LET  F=SOP  (S/PE) 

509  SCPRTCH  *3 

510  GOSUB  700 

512  IF  V5=  0 GO  TO  56  0 
514  PPIMT 
516  PPIMT 
520  PPIMT 

522  PPIMT  * M ="M 

524  PPIMT  ' H ="H 

525  PPIMT  " DEL. CM  ="B<H)*2.54 

530  PPIMT  " TRUW.PSF  =“S 

535  PPIMT  • UTRU.M/S  ="F*.3048 

536  PPIMT  " U<D).H/S  =“U9/3.2808 

540  PPIMT  " PI  =" <l/4.86)*< <U4/F)-2.43*L0G<R3)-5) 

550  PPIMT  “ PMS  =“Z  <M> 

551  LET  M2=M1 — 1 

552  LET  Ml=l 

553  PPIMT 

554  PPIMT 

555  PPIMT 

556  GOSUB  700 

557  GO  TO  600 
560  LET  G<H>*Z<'M> 

562  IF  H=>£  GO  TO  570 
564  LET  H=H-M 

566  LET  B<H)*B<H-1)+B1*>B<H-1> 

568  GO  TO  285 

570  IF  G«H)>G<H— 1)  GO  TO  576 

572  IF  B <H) >B  <H-1)  GO  TO  586 

574  GO  TO  580 

576  IF  H=2  GO  TO  580 

578  GO  TO  592 

580  LET  H-H+l 

592  LET  B<H)*B<H-1)-B1*B<H-1> 

584  GO  TO  285 
586  LET  H=fKl 

588  LET  B < H) =B  <H— 1> +B1 ♦B  <H— 1 ) 

590  GO  TO  285 

592  LET  B <H>  =B  <H— 1> 

593  LET  V5=l 
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5^4  60  TO  490 
600  LET  Nl=*8 

610  LET  N2=N3 

611  PRINT 

612  LET  V6= 1 

620  GOSUB  700 
624  LET  R2=U5<N> 

626  LET  U7=U6<N> 

630  FOP  N=1  TO  N3 

632  LET  D4=D4+  <1  - > <U  <N>  -HJ  <H+ 1 > ) ♦.  5/U9>  ) ♦<  Y <N>  -Y  <N-1  > >*2.54 

634  LET  D5=R2-<  <U5  < N> -HJ5  < N-l  > > ♦.  5> 

636  LET  D6=D6-*-D5*<  Y <N>  — V <N- 1 > > *2.  54 
638  LET  D8=D8+D5*D5*-<  Y (N'*  -Y  •'N—  1 > ) ♦£ . 54 
64  0 IF  B>H»<=Y<'N>  60  TO  644 
642  NEXT  N 
644  PRINT 

646  PPINT 

647  PPINT 

648  PPINT  " INTEPGRL  PPOPEPTIES  OF  BOUNDRPY  LRYEP" 

649  PPINT 

650  PPINT  ‘ DELTR  STRP  K <CM>="D4 

652  PPINT  * DEL  •:CN'>="D6 

654  PPINT  " 6 =“D8^D6 

656  PPINT 

658  PPINT 

660  PPINT 

661  SCPRTCH  =2 

662  PPINT  ■ VELOCITY  DEFICIT  COPPELRTION- 

664  PPINT 

666  PPINT  -Y/DEL”.  “U*DEF  CRL-."U*DEF  NERS" 

668  PPINT 

670  FOP  N=1  TO  N3 

672  IF  B<H><=  Y<W  GO  TO  676 

674  PPINT  Y<N4€.54/D6»U7-U6(N> »R2-U5(N> 

675  WRITE  s2»Y^N>*2.54/D6»U7-U6<N> *R2-U5<N> 

676  NEXT  N 
678  STOP 
700  LET  T8=0 

705  FOP  N=1  TO  N3 

710  IF  B <H>  > Y (N>  GO  TO  74  0 

72 0 LET  IJ4=U3  <N- 1 > + <U3  <N>  HJ3  <N- 1 > ) *<B  <H>  -Y  <TN- 1))/  <Y  -Y  (N- 1 > > 
722  LET  U9=U  < N- 1 > ■*-  Ol  <N>  -U  <N- 1 > ) *<B  < H)  -Y  <N- 1 > ) ✓ <Y  <N>  -Y  < N- 1 > > 

730  GO  TO  750 
740  NEXT  N 

750  LET  R3=B  <H) *F/ < 12*V> 

760  FOP  N=N1  TO  N2 

770  LET  W=2*<'  <S  IN  < a . 5707*Y  <"N>  > / <B  <W  > > > ^2> 

780  LET  L=. 5*<  CU4/F)  — <2. 43*L06  <R3) ) —5) 

790  LET  L 1=2. 43*L06  <"Y  ^N>  ♦F/'  < 12*V>  ) 

800  LET  U1«L145+L*W 
810  LET  U2=U3  <N>  sf 
851  LET  Wl  = <U2-l_l-5>  /L 
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852  LET  U5 (M>  *02 

853  LET  U6<M>»U1 

860  LET  T8='U1-U2:>~2«-T8 
870  IF  V5  =0  GO  TO  893 
872  IF  B<HX*Y<M>  GO  TO  910 

880  PPIMT  Y (M>  »2 . 54 » Y <MXB  <H>  , Y <M!>  ( 1 2^V>  *U1  ,U2 

881  IF  V6=l  GO  TO  885 

882  IF  Ml  = l GO  TO  885 

883  PPIMT  TRB  <15> * W1 

884  PPIMT 

885  IF  M*>K9  GO  TO  892 

886  LET  Wl=3 

887  IF  <Y(M)«f/(12»V))*>  10  GO  TO  892 

888  LET  Ul*<'Y<M>*F/'<ie*V>> 

892  WPITE  “3  * Y (M>  XI  2*V>  * U 1 * U2  * Y (M')  ✓ B <H>  » W 1 

893  MEXT  M 

900  LET  Z (M>  =SQP  <T8) 

910  PETUPM 

1000  BRTR  "3MFL-125"* ‘3VDE-125-* "3VCO-125" 

1010  DRTR  -1.25*21. 8. 540*. 68*. 75. 24, 69*100. .01,. 01 
9999  EMB 


RE«rr 


PROGRAM  TBLJNDIM 


010  PEM  TBLJNDIM 

11  PEM  P9  IS  PRMP  NO.Z  IS  NO  OF  X STA.  X<L> 
1£  PEM  DELTA  > CM ^ FPOM  VCOLES.  P5<L>  IS  PWALL 
100  FILES  F15F2 

110  DIM  Y'100>  .U'100>  .Raoo>  .P<100>  »TO<100> 
120  DIM  PI  '.£>  *U1  < £> 

130  DIM  X > 1 n . D < 1 1 ) » P5  (1 1 > 

140  PERD  P^.Z 

150  FOP  L= 1 TO  11 

160  PERD  X <L> * P <L>  * P5  <L> 

170  LET  P5 -L> =P5 'L> *97309 
190  NEXT  L 

190  PRINT  “ MERN  FLOW  PPt 

£00  PRINT  * PRMP  N0'P9 


IS  XSTA* IN>  . D <L) 
iPW^PO) 


MERN  FLOW  PROFILES' 

PRMP  N0'P9 

X STRT  I ON  * CM=  “X  '.Z^  ♦£  .54 


F 1 S « F £S 


£10  PRINT  " X STATION. CM=“X'Z>*£. 54 

££ 0 PRINT 
£30  PRINT 
£40  PRINT 

£50  PRINT  * WHAT  RPE  FILES  NEWFLOWX . MERNXXXX' : 

£6  0 INPUT  FIS*  F£S 

£6£  PPINT 

£64  PPINT 

£66  PPINT 

£70  FILE  si. FIS 

£80  FILE  =£.F£S 

£90  : CRATCH  s£ 

300  FOP  L= 1 TO  Z 
310  FOP  M= 1 TO  100 
3£ 0 IF  L=Z  GO  TO  350 
330  PERD  =1 .&1 .G£.G3.Q4.Q5 
340  GO  TO  360 

350  READ  si  »Y  (M>  .U'M'1  » R <M>  »P<M)  » T 0 <M"> 

35 £ LET  P ■ M"1  =P  f M>  *1 33.  3 

360  NEXT  M 

370  NEXT  L 

400  PPINT 

410  PPINT 

4£ 0 LET  P ■ 0>  =.  00349*95  <Z~>  /300 

430  LET  M4=SGP • 5*< < <97309'P5 > Z> > £857> -1> > 

44  0 LET  T 4=3  1 7 d *.  £*M4*M4’> 

450  LET  P4=.  00349*P5<'Z)^T4 
460  LET  IJ4=M4*£0.  04*SGPa4> 

470  LET  PI  < 1>  =P4 

49  0 LET  U 1 1 1 > =1.14 

490  FOP  N= 1 TO  100 

500  IF  Y (N^  > D CZ>  GO  TO  5£0 

510  NEXT  N 

5£0  LET  U£=U  > N-H  * ('U  <N>  -U  <'N-1)  > *<  D <2^  -Y  <N-1>  > ✓ <Y  (N'1  -Y  <N-1>  > 

530  LET  U£=U£/'.  995 

540  FOP  N= 1 TO  100 

550  IF  U'N>>U£  GO  TO  570 

560  NEXT  N 


46 0 LET 
470  LET 
490  LET 
490  FOP 
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LET  Y2=Y  CJN—  1 ) 4-  ^ Y ( N'*  -Y  -M-  1 > > ♦ < U2-LI  <M-1>  > / <U  1 M'1  —LI 1 M—  1 > ) 
LET  Y3=  < Y2-Y  <M-1>  > / <Y  W -Y  <M-1>  ) 

LET  P2=P<M-1>+<P'M>-P  <M-1>>^Y3 
LET  'M-1))»Y3 

LET  T3=T 0 <M—  1 ) ♦ <T0  <M>  -TO<M-D)  ♦YS 
LET  93=20.  04^S9P  <T3> 

LET  T2*T3^<1-.2^<U2*02/'  <93«*3>  > > 

LET  M2=SGP  <?♦■<  CT3/T2>  — 1)  ) 

LET  V2=.  OOOOOl  0969^-<  <T2-*-l . 8)  . 5>  / < <T2M . 8)  +199) 

LET  S2=LI2«4?2/V2 

PPIMT  " YE. CM  ="Y2 

PPIMT  “ TOE. DEG  K ="T3 

PPIMT  - TE.DEG  K = "T2 

PPIMT  “ PE  * M ''M2  =“P2 

PPIMT  " IJE.M/S  = "U2 

PPIMT  - PHOE.KG/M3  = "P2 

PPIMT  - PEE.M-1  = "S2 

PPIMT  - ME  = "M2 


U''UE  " . “PHO/PHOE- 
. "TO-'TOE".  “P-'PE’ 


LET  V2=.  000001  0969^<  < T2^1  . 8>  ~ 1 . 5)  / < <T2^1 . 8>  +1  9 
LET  S2=U2«*?2/V2 

PPIMT  ' YE. CM  ="Y2 

PPIMT  “ TOE. DEG  K ="T3 

PPIMT  - TE.DEG  K = "T2 

PPIMT  “ PE.M/M2  ="P2 

PPIMT  " UE.M/S  = "U2 

PPIMT  - PHOE.KG/M3  ="P2 

PPIMT  - PEE.M-1  = "S2 

PPIMT  - ME  = "M2 

PPIMT 
PPIMT 
PPIMT 

PPIMT  “Y"  • " Y''D " . “M _ . “U/'UE”.  “PHO/PHOE" 

PPIMT  "PE'  PEE"  . "T/TE"  . "TO/TOE".  "P-'PE".  "H" 

PPIMT 

PPINT 

PPIMT 

FOP  M=1  TO  100 

LET  9=2  0.  04»SQP  'TO  ffO  > 

LET  T=T  0 <M>  +■<  1 - . 2 ♦Ll  <M>  ♦Li  <M)  ✓ <9^9\>  > 

LET  M=C0P<5^'  <T0<M>xT>-l>> 

LET  V=.  0 000 0 1 0969^ < <l+-l . 8>  ~ 1 . 5>  / < <TM . 8>  ♦ 1 99> 

LET  S=P  <M^ ♦U  <M> /V 

LET  H=<T0>M>-3  00>/ CT3-300) 

PPIMT  Y <M>  . Y <M>  /D  <Z>  • M.  U <M> /U2.R  <M>  ^R2 
PPIMT  S/S2. T/T2. TO  <M) /T3.P  <M> /P2.H 
PPIMT 

WRITE  -2. Y <M> /D  <Z> . M*  U <M> /U2. P <M> /P2. H 
IF  Y <M> > Y2  GO  TO  930 
LET  W=Y  <M>  — Y <f+—  1 > 

LET  Ml =SQP  C5+<  < <973 09/P  <M>  > ^ . 2857> -1 > > 

LET  Tl=31 7/ < 1 ♦. 2^M1 *M1 ) 

LET  PI  <2)=.  00349^<P<M>/T1> 

LET  U 1 <2>  =M  1 +2  0 . 04*SG>P  <T  I ) 

LET  Cl=Cl  + < <U1  <2>-MJl  <l))/2-<U(M)-H.KN-l))/2)^ 
LET  V=  <P  <N>  ^Ll  <M>  ♦P  <M—  1 > ♦€!  <M—  1 ) ) /2 
LET  V1  = <P1  <2)  ♦Lil  <2>  ♦P 1 < 1 > *0 1 <l>>/2 
LET  C2=C2-*-  <V  1 — V>  ♦W 

LET  V=  <P  <N>  *0  fM>  ♦L*  <N'>  ♦P  <H—  1 > ♦Ll  <M-1  > ♦LI  <M-1  > > ✓2 
LET  Vl  = <P  <"N>  ♦Li  <M>  ^Lll  < 2)  -H?  <M— 1 ) «*!  <M- 1 > ♦Ul  Cl))^ 
LET  C3=C3+<Vl-V>«*l 
LET  PI  <1>=P1  <2) 

LET  Ul  <1)=<J1  <2) 

MEXT  M 
PPIMT 


<*♦*>> 
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042  PPIMT 
Q44  PPIHT 
*46  PPINT 

*49  PPIHT  “ DELTR.CM  ="D<Z> 

950  PPIHT  " DEL«+ -Cn  ="C1-U4 

953  PPIHT  “ DEL*. C(i  = "C2-  ' U4*P4  > 

•^54  PPIHT  " THETR.Cn  = "C3-  UJ4*U4*P4T 

956  PPIHT  “ PE  THETR  ="32*C3  • 1 0 0*d4*Li4*P4) 

*70  STOP 

1000  DRTR  3.2 

1010  DRTR  0*  . 7042, . 02986 

1020  DRTR  .25, .6835. . 03212 

1 03 n DRTR  . 75, . *566, . 04 

1040  DRTR  1 .25*  . b767.  . 04732 

1050  DRTR  1 .5. .6632, . 05361 

1060  DhT R 1 . 75. .6632. . 05955 

1070  DRTR  2. .6835, . 06504 

1080  DRTR  2.25, .6903, . 06899 

1090  DRTR  2.5. .6835, . 07207 

1100  DRTR  2.75* .71 12, . 07882 

1110  DRTR  3, .7401, . 08277 

PEftDY 


I 
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PROGRAM  RAMSHER 


10  PEM  PBMSHEP 1 

11  PEM  H IS  PBMP  MO.  N 1 IS  NO  OF  FIRST  X STB  USED  IB  CURVE  FIT.  N2  IS 

12  PEM  THE  NO  OF  THE  LBST  X STB  IN  CUPVE  FIT,  Z IS  THE  X STB  FOR  TBU  CBLC 

13  PEM  ’MUST  BE  <N2> • W IS  PEFEREHCE  X STB  WHICH  DETERMINES  Y VBLUES. 

14  PEM  K 1 IS  NO  OF  TIMES  1ST  X STB  INCLUDED  IN  CUPVE  FIT  <USUBLLY*1> . 

15  PEM  F <N'>  IS  TBUW-MMHG> .R5<N>  IS  RBD  OF  CURVBTUPE  <XM>  » X <N>  IS  X STB « IN 

16  PEM  Bl<N->  IS  SUPFBCE  SLOPE  <RBD> , D6<N>  IS  LOCBL  HEIGHT.  Y <CM>  » OF 

17  PEM  SUPFBCE  WITH  RESPECT  TO  PBMP  L.E. 

100  FILES  F15F2.F3 

110  DIM  YO  1 » 1 00>  ,U<  1 1 » 1 00>  .P<  1 1 . 1 00  ,P<1 1 . 1 00) 

120  DIM  Cl  <1 1>  »C2<1 1)  ,C3<1 1)  ,C4(1 1>  ,C5(11)  ,S  <1 1> 

130  DIM  B <3*  3)  , B < 3, 1 > » D1 <2>  » D2  <2>  » D3  <2>  » D4  <£> ,X1 <1 ,3> 

140  DIM  P5  < 1 1)  ,F  <1 1>  »X  ai>  »B1  <1 1>  ,D6  Cl  1>  . U1  <3»3>  »Z6(3»  1> 

142  DIM  C6 1 1 1 > 

230  PEBP  H.N1 . N2.Z. W.Kl 
240  FOP  N=1  TO  N2 

25  0 PEBD  F ^N> . R5  <N>  » X <N>  » B1 <N>  » D6  <N> 

252  LET  X <N> =X  < N> *-2. 54 
254  LET  F <N>  *F  <N>  ♦-!  33. 29 


260 
270 
28  0 
290 
292 
294 
296 
300 
31  0 
312 
314 
316 
320 

322 

323 

324 

325 
330 
34  0 
350 
352 
36  0 
370 
380 
382 

384 

385 

386 

387 

388 

389 

390 


TUPBULENT 


X 


WHBT  BRE 
F1S.F2S.F3S 


FILES  NEWFLOWX. SHERXXXX » EDYXXXX" 


NEXT  N 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
PRINT 
I NPUT 
PRINT 
PRINT 
PRINT 

FILE  si. FI* 

FILE  s£,F2* 

FILE  S3.F3S 
SCRBTCH  ~2 
SCPBTCH  *3 
FOP  N«N1  TO  N2 
FOP  M=1  TO  100 

PEBD  sl.Y(N.M) »U(N,H>  »P(N»M> .P<N,M> ,09 
LET  P <N»  M>  *P  <N»  M>  »133. 29 
NEXT  M 
NEXT  N 

PRINT  -YCCMl  " , “ TBU/TBLtW"  .“I1"»“I2“»"I3 

PRINT  TBB<15>, “14- 

PRINT 

PRINT 

LET  D1 <l>-0 
LET  D2<n-0 
LET  D3<1)*0 
LET  D4  < 1 ) *0 
FOR  M«1  TO  100 


SHE BP  STRESS 
PBMP  NO“H 
STBTIDN <CM>-“X  <Z> 
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4 0 0 
410 
42  0 
4 3 0 

44  0 

45  n 

46  0 
470 
472 
474 
476 
478 
48  0 

482 

483 

484 
48* 
488 

4 9 0 
482 
4^4 
4^5 
4C** 
497 
4Q8 

5 n fi 
5 0£ 

510 
515 
52  0 
525 
530 
535 
54  0 
545 
550 


56  ii 
565 
567 
570 
578 
58  0 
581 

585 

586 
58  7 
5Q  n 
5QK' 
6 on 
* 08 
*05 


LET  Y*Y<W,H'» 

LET  00001 

FOP  L= 1 TO  N2 

LET  Z2=Y*C0S<81  *L>>-*-B6(L> 

LET  Z1=Y*C0S<81  <L-n>-*-P6tL-l> 

LET  L2=X>L'*-Y*7IN>81  iL"*  > 

LET  L 1 =X  <.L- 1 > — Y*S  I N < R 1 « L — 1 > > 

LET  23=72-71 
LET  l3=L2-L 1 

LET  7 'L^sSOR"  vL3*L3''  ♦ ^23»73''  > +S  <L-1> 

ME  XT  L 
M8T  R=ZEP 
LET  * = 0 

FOP  L=N1  TO  M2 
LET  ► =*>  -M 

LET  8-  1.  P=M<  1.  1>-M 

LET  8 <■  1 • 2>  =H  ' 1 ? 2> ♦S  <L>  CL-* 

LET  8 ' 1 • 3^  =8  > 1 • 3 ♦S'  <L  > 

LET  8 > 2 • 2>  =8  • 2*  2>  ♦ £ <L>  ^L')  ♦S  •■L'*  ♦£  «.L) 

LET  8 ' 2«  S'*  =8  1 2*  S'*  ♦ $ < L '•  • L"*  *7  <L> 

LET  8 • 3 • 3'*  =8  < 3 - 3:*  ♦ $ <L>  iL> 

IF  L>1  60  TO  4^7 
IF  K < ► 1 60  TO  483 
MEXT  L 

LET  8 '2. 1) =8< 1 »2) 

LET  8-3. 1>=8<1*3> 

LET  8 <3, 21)  =8  <2.  S'* 

FOP  L =M 1 TO  N2 
FOP  M=1  TO  100 
IF  Y <L « M*  ='*  Y 6D  TD  525 
NEXT  M 

LET  Y 1 = - Y-Y  • L • M-l  ^ •:  Y <L  . M~*  -Y  <L  * N-l  > > 

LET  Cl  >L  * =U  *L«M-n  ♦ »L -M^ -U  > L»M-1>  ■*  ♦Yl 
LET  P 1 =U  * L * M- 1 > ♦P  cL» M— 1 ) 

LET  F2=U  * L • M'*  *P  • L * M' 

LET  C2  <L>  =F1 ♦ <B2-P1> ♦Yl 

LET  F 1 =U  1 L • M- 1 > *.L * M- 1 > <L  * M- 1 > 

LET  B2=LUL.N>*Li'L.N'**F*<L»N> 

LET  C3*L'*  =B  1 ♦ * F2 - B 1 > ♦>'  1 

LET  C 4 • L > =P  * L * M-  1 > ♦ •:  P ( L » N ^ — P • L • M- 1 ) > ♦Y 1 
LET  C 6 * L~>  =P  * L ♦ M-  1 > ♦ > P <.  L • N'*  — P < L » M-  O'*  ♦Y  I 
NEXT  L 
LET  K=0 

FOP  L =M 1 TO  N2 

LET  ► =R  ♦ 1 

LET  C5  ■ L'*  =C1  <L> 

IF  L> 1 60  TO  5Q0 
IF  ► <>  1 60  TO  581 
NEXT  L 
60SUB  800 

LET  PI  *2'*  =X  1 * 1 . n+Xl  < 1 .3^  »S  c2>  ♦Xl  < 1 »2>  ♦S  • Z'**?:  'Z> 
LET  *>=o 

FOP  L =N 1 TO  N2 


I ss 


606 
61  0 
611 
612 
615 
62  0 
625 
629 

63  0 
63 1 

635 

636 

637 

64  0 
645 
650 

654 

655 

656 
66  0 
661 
662 
665 
670 

672 

673 

674 

675 

676 

677 
673 
68  0 
685 
690 
695 
700 
705 
707 
71  0 
715 
720 
725 

73  0 
735 

74  0 
"’45 

747 

748 

749 

750 
752 


LET  K=**l 
LET  C5<L)=C2<L> 

IF  L> 1 GO  TO  615 
IF  ¥ <¥ 1 GO  TO  606 
NEXT  L 
GOSUB  800 

LET  B2  <2>  =X  1(1. 3>  «-2*X  1 <1 . 2>  <Z> 

LET  K=0 

FOP  L=N1  TO  N2 
LET  K=K-*-l 
LET  C5  <L>  =C3  (L) 

IF  L> 1 GO  TO  640 
IF  ¥<¥1  60  TO  631 
NEXT  L 
GOSUB  800 

LET  D3  <2>  =X  1 •:  1 , 3>  •*-2*X  1 <1  * 2>  *-S  <Z> 

LET  K = 0 

FOP  L=N 1 TO  N2 

LET  ►'=♦<>  1 

LET  C5a>=C4CL) 

IF  L> 1 GO  TO  665 
IF  ¥ <¥■ 1 GO  TO  656 
NEXT  L 
GOSUB  800 

LET  D4  <2>  =X  1 a » 3>  ♦ 2*X  1 < 1 » 2)  *S  <Z> 

FOP  L=N1  TD  N2 
LET  C5 <L> =C6 <L> 

NEXT  L 
GOSUB  800 

LET  D5=X  1 a.  1>+X1  < 1 » 3>  ♦£  <Z>  ■♦■XI  <1  , £>  ♦£  <Z)  *S  <Z> 
LET  K=l/P5  <Z> 

LET  T=Y.:W.H>-Y<W,H— 1> 

LET  Bl  = l/  a «+*•<  <Y<W,M>-*-Y  <W,M-1> > /2> > 

LET  D=  <D3  <2>  ■♦'D3  (1)  )/£ 

LET  11  = 11  -*-B  1 ♦Ii^T 

LET  D=(D2<2>-*-D£<l)>/e 

LET  V=  <D1  <2>  4-D1  (1>  ) ^2 

LET  G1=<D1 <2>-Bl <!>>/T 

LET  G=G-*-B*T 

LET  I2=V*B1*G 

LET  I3=I3-*-2*I2*Bl«**T 

LET  B='B4<2)+ri4a>>/2 

LET  I 4= I 4+B 1 ♦B^T 

LET  Tl  = H-<l/F<Z>>*<Il-I2-I3-*-I4> 

PPINT  Y • T 1 > Il^F-'Z)  » I2/F<Z)  »I3/F(2D 
PPINT  TRB<15) * I4/F CZ) 

PPINT 

WPITE  ~2fY*Tl»D£<2) »D3(2) »B4<2) 

WPITE  “3 » Y»  T 1 *F <Z> » G1 * D5 
LET  Dia>*ri<2) 

LET  D2 (1>  =D2  <2> 
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754  LET  r3^n=r3‘£'> 

75*  LET  P4<1>=P4<3> 

760  NEXT  M 

765  7TDP 

800  NRT  F=ZEP 

805  EDP  L =N 1 TD  N3 

810  LET  F< 1. 1>=B<1» 1>+C5<L> 

815  LET  F(3,  1>=F<3.  1 > +S  <L>  <L>  «C5  <L> 

830  LET  F <3«  1 > =F  c3»  1)  +S  <L>  *C5  <L> 

835  NEXT  L 

830  MRT  W1=CDN<3»3> 

835  MRT  W1=1NV<R> 

840  MRT  Xl*CnN<l»3> 

845  MRT  Z6=CDN<3,1> 

850  MRT  Z6=W1*B 
855  MRT  XI = TPN  <Z6> 

860  PETUPN 

1 000  PRTR  3.4.11 .6.7.  1 
1010  PRTR  ,36.-1000000.  0.  0.0 
1030  PRTR  .363* -18. 14. ,£5. . 03748. . 01364 
1030  PRTR  ,376. -33. 5« . 75. . 1 005* . 1 036 
1040  PRTR  .346.-36. 75. 1 .35. . 1537. .3659 
1050  PRTR  .307.-38.84, 1.5. . 1759. .3716 
lOtrvO  PRTR  . 336.— 30.  91.1. 75.  . 1976.  . 4919 
1070  PRTR  .347.-33.97.3, .3179. .6361 
1080  PRTR  .376,-35. 03-3.35. .3371, .7734 
1090  PRTR  . 393. -37. 07. 3. 5. . 8553. . 9333 
1100  PRTR  . 486,-39. 1 1 ,3. 75. . 8736, 1 . 1 053 
1110  PRTR  . 44^. -41 . 55, 3. . 3891 , 1 .3886 

QQQCl  ehp 
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PROGRAM  EDDY 


10 
1 1 
12 

13 

14 
1 0 0 
1 1 0 
190 
200 
210 
22  0 
230 
232 
234 
236 
238 

24  0 

25  0 
260 
270 
28  0 
282 
284 
286 
290 

295 

296 
300 

31  0 
312 

32  0 
330 
340 
350 
352 
354 
356 
360 
400 
1000 
10  07 
1008 
1 009 
1010 
1 020 
1 030 
1 040 
1050 
1060 
1070 
1 080 
9999 


STR.X<L)  IS 
SUP  K •D**'’ 
UE  vR/SEC) 


RRE  FILES  EDYXXXX * V I SCXXXX “ • 


PER  EDDY 

PEn  P8  IS  THE  PRRP  ND.Z  IS  THE  HD  DF  THE  X 
PEH  D-L)  IS  DELT8  <CH).D1<L)  IS  DELTR  STRP 
PER  METERS  FRDM  VCDLES  DP  TFLJNDIR.  LKL)  I 
PER  TBLJHDIM 
FILES  FI J F2 

DIM  X <.  1 1 ) * D < 1 1)  , D1  U 1 ) . LU 1 1) 

PERD  P8»  Z 
FDP  L= 1 TD  Z 

PERD  X <L>  » D (L)  * PI  CL)  » U CL) 

HEXT  L 

PPIMT  ’ EDDY  VISCDSITY  CRLCULRTIDH 

PPIMT  “ PRMP  MD"P8 

PPIMT  “ X STRTION="X (Z) ♦£. 54 

PPIMT 

PPIMT 

PPIMT 

PPIMT  “ WHRT 

I MPUT  FIS.  F2$ 

FILE  ::1.F1S 
FILE  S2.F2S 
PPIMT 
PPIMT 
PPIMT 

SCPRTCH  -2 

PPIMT  " Y “ • “ Y/D’ 

PPIMT 

FDP  M= 1 TD  100 
PERD  “l.Y.T.G.P 
IF  T < 0 GD  TD  354 
LET  L=SOP  <T /P) / <D  <Z) ♦€) 

LET  V=T / <U  <Z)  *B  1 <Z)  ♦RoG) 

PPIMT  Y » Y/D  <Z>  * L*  V. G^D  <Z)  /'Ll  <Z) 

“2 * Y.  Y/P  <Z)  »L.V,G*TKZ)/U<Z) 

36  0 

Y.Y/DCZ)  . TRP  '61 ) !G-*D<Z)/LKZ) 

~2. Y. Y/D  <Z) • 9999. 9999. G^D <Z> /U <Z) 


L'D“  ■ 


E/UEP^K" ? " D 1 U - 


WRITE 
GD  TD 
PPIMT 
WRITE 
NEXT 
STOP 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
DRTR 
EMP 


M 


3.4 

0.  1.  1 . 1 
.25.1.1.1 
1.  1 . 1.  1 

1.25. . 6767.. 0944.610.2 

1.5. . 6632. . 1 19.605 

1.75. . 6632. .112.600.5 
2.-6835. . 1138.596.4 

2.25. . 6903.. 1069.587.6 

2.5. . 6835. .1021.581.4 

2. 75..  71 12.. 1045.580 

3. . 7401. . 1033.573 


X STR> IN’ 
IN  CENT  I - 
FROM 


UE)  / D Y/D) 
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APPENDIX  B 

COMMENTS  ON  THE  MIXING  LENGTH  CONSTANT 

In  Section  4.6,  it  was  shown  that  in  the  region  of  adverse  pressure  gradient 
(APG),  the  constant  x in  the  mixing  length  relation  L = xy  was  0.65  in  con- 
trast to  0.4  for  the  zero  pressure  gradient  (ZPG)  flow.  This  implies,  there- 
fore, an  inconsistency  of  the  APG  data  with  the  Law-of-the-Wall  correlation, 
since  the  latter  is  based  on  the  value  k * 0.41,  and  raises  the  question 
whether  the  increase  in  x is  real  or  possibly  due  to  experimental  error. 

Recall  that  the  derivation  of  the  wall-wake  velocity  correlation 

u ■=  - tn  y+  + C + — W (y/6)  B.l 

X H 

includes  the  assumption  l = xy.  References  8 and  11  document  the  basis  for 
selecting  the  currently  accepted  values  of  the  constants  ;i  = 0.41  and  C = 

5.0.  These  values  are  concluded  to  be  independent  of  pressure  gradient, 
whose  effects  are  reflected  instead  in  the  magnitude  of  the  wake  strength 
parameter  “.  It  should  be  noted  that  the  values  of  x and  C are  based  on  data 
for  which  100  < y+  < 300,  since  it  was  felt  that  closer  to  the  wall  experiment- 
al errors  may  cause  the  measured  velocity  to  be  too  high,  while  farther  away 
the  effects  of  the  wake- like  outer  flow  become  large. 

In  the  present  study,  the  curve-fit  described  in  Section  4.2  was  restricted 
to  data  for  y >50.  A plot  of  the  velocity  profiles  in  u+  - y+  coordinates, 
shown  in  Figure  28  for  Ramp  3,  is  reproduced  in  Figure  B.l  which  includes, 
now,  lines  of  constant  y/6.  Notice  that  for  the  ZPG  case,  y/6  = 0.05  corresp- 
onds to  y+  ~ 20.  Moving  downstream  in  the  APG  region  the  v+  value  correspond- 
ing to  y/6  = 0.05  increases  (the  increase  in  Reynolds  number  shifts  the 
velocity  profile  to  the  right,  i.e.,  to  larger  y+  values)  so  that  at  the 
rear  of  the  ramp  y+  at  y/6  = 0.05  is  about  50. 

A plot  of  1/ 6 versus  y/6  to  an  enlarged  scale  is  shown  in  Figure  B.2  which 
indicates  clearly  chat  the  slope  x shifts  rapidly  from  0.4  to  0.65  along  the 
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ramp  and  that  the  linear  portion  of  the  1 versus  v variation  extends  from 
y/6  *=  0.05  to  almost  0.2  (see  also  Figures  63  and  64  to  define  the  upper 
limit).  A value  of  y/5  = 0.05  corresponds  to  about  0.035  cm  which  is  near 
the  outer  edge  of  the  viscous  sublayer  while  v/5  of  0.2  corresponds  to  y+ 
ranging  from  80  to  200  for  the  profiles  in  Figure  B.l.  Therefore,  most  of  the 
data  points  (particularly  for  the  APG  region)  shown  in  Figure  B.2  are  included 
in  the  curve-fit  to  the  Law-of- the-Wake , Equation  4.5.  However,  these  points 
are  small  in  number  compared  to  those  in  the  outer  portion  of  the  boundary 
layer,  where  contributions  from  the  wake  function  are  large,  and  although 
they  deviate  from  Equation  4.5  (i.e.,  they  reflect  = 0.65  instead  of  0.4) 
the  rms  error  of  the  curve-fit  is  still  small. 

For  both  ramps,  the  u+  - v+  profiles  for  representative  stations  located 
upstream  of  the  leading  edge  (ZFG),at  its  mid-point,  and  at  the  rear  station 
have  been  plotted  in  Figure  B.3  together  with  the  sub-layer  relation  u+  = v+, 
Equation  4.8,  and  the  Law-of- the-Wal 1 , Equation  4.4.  It  is  observed  that 

1)  For  the  2PG  case,  the  data  for  y+  < 50  fall  on  or  parallel  to  Equation  4.4 
until  v+  20  (v/5  = 0.05)  where  they  then  merge  with  the  sublayer  profile. 

2)  For  the  APG  region,  the  data  below  y+  = 100  deviate  from  Equation  4.4, 
showing  a smaller  slope  (i.e.,  a larger  value  of  ;t)  than  the  Law-of- the-Wal 1 
until  v+  ~~  25  where  again,  they  approach  the  sublayer  profile.  The  trend  of 
the  APG  data  is  clear  and  consistent.  In  fact,  since  = 0.65  is  determined 
directly  from  the  mixing  length  calculations,  this  guarantees  that  the  APG 
data  will  not  fit  the  classic  Law-of- the-Wal 1 (with  . = 0.41). 

In  view  of  this  discussion,  the  possibility  of  experimental  error  causing  the 
observed  increase  in  in  the  APG  region  is  considered  unlikely  since: 

1.  The  effect  is  not  observed  in  the  ZPG  case  although  the  techniques  and 
instrumentation  were  the  same  at  all  x stations. 


I 


1 b0 


2.  The  data  tor  all  stations  w«r« 

boundary  layer,  that  is,  tor  \ | 

the  wall  tor  probe  errors  to  . • . 

3.  Sturek  obtained  the  s r 

9^.  Thus,  two  experiments  * mu!  . 

operating  conditions,  wit: 
yield  the  same  result. 

It  appears,  then,  that  th«  . 

APG  flow  is  not  due  to  c\p« i i!*  • i 

from  the  ZPC.  value  is  not  Know 
results.  It  is  suggested  t • . 
gradient  or  longitudinal  lurv.iiun  i,, 

there  is  a need  tor  more  inform. 1 1 i or 

It  is  instructive  here  to  apple  Spaldin.  ' • > « - 

in  the  buffer  region  to  the  APG  data  m otdet  • a >1.  . . . . - 

velocity  profile  associated  with  • 0 . t-  ■ . Tit  j - . not  a - y i.  t lor 

4.8  to  the  point  y+  * 100,  u+  “ lb. 2 where  the  data  fro;-  -ot  i >\  rid 
APG  profiles  agree  and  are  identical  to  the  conventional  . -tl-the-k.  11, 
Equation  4.4.  At  this  point,  we  assume  *=  0.b5  and  use  Iquatiot  ...8  to 
calculate  K-, . With  both  and  known,  we  can  flier,  use  Equation  -*.8  to 
calculate  u+  versus  y+.  The  results  are  plotted  in  Kigurt  P . J where  they 
seem  to  coincide  with  the  APG  data  throughout  the  entire  y+  range  near  the 
wall . 
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